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An investigation has been made of certain phenomena 
associated with space charge external to the grid which 
surrounds the cathode in a multi-grid tube. This paper 
reports a new effect, the essential feature of which is a 
negative plate current-emission characteristic. The con- 
ditions under which this phenomenon is observed are out- 
lined. Results of experiments with double and triple grid 


tubes are summarized. In general, two points on the 
characteristic may be found for which the plate current is 
independent of small variations in emission, whether of a 
statistical or “‘functional’’ nature. An explanation is 
presented based on the interaction of multiple space charge. 
An application to the problem of high gain amplification is 
suggested. 





I. INTRODUCTORY 


HE study of the effect of space charge on 

thermionic currents has been restricted 
mainly to those cases in which maximum density 
of electrons occurs in close proximity to the 
emitter. Under certain conditions, control of the 
current in multi-grid tubes may be produced by 
what we shall designate as‘‘ external” or ‘‘mul- 
tiple’ space charge. These terms are applied 
here to space charge resulting from accumulation 
of electrons in regions of the tube external to the 
grid which surrounds the cathode. One case in 
which this situation was used has been treated 
briefly by the writer in connection with studies 
of the shot effect.!: ? 

The characteristics associated with external 
and multiple space charge have been investi- 
gated. It is the purpose of this paper to report a 
new effect which appears under these conditions. 


II. GENERAL DISCUSSION 


Space charge limitation in the conventional 
sense involves a saturation level beyond which 


1E. W. Thatcher and N. H. Williams, Phys. Rev. 39, 
474-496 (1932). 
?E. W. Thatcher, Phys. Rev. 40, 114 (1932). 
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the current ceases to rise with emission. The 
essential feature of the effect reported here is a 
situation in which the plate current decreases 
with increasing emission, preceded and followed 
by regions in which the current rises with 
emission. The curve which represents the plate 
current as a function of emission will have two 
horizontal tangents and a region of negative 
slope. 

Experiments with double and triple grid tubes 
have shown that the conditions favorable to the 
appearance of this effect include : 

(1) An accelerating field for electrons as they 
leave the emitting surface. This is maintained by 
holding the first grid at sufficiently high positive 
potential to prevent space charge limitation of the 
current from the cathode. This also serves to 
prevent the escape of positive ions, which are 
usually, if not always, present in thermionic 
emission.’ 

(2) A retarding field in which the stream 
velocity is reduced before passing into the outer 
space. This is produced by making the second 
grid only slightly positive with respect to the 
cathode. It may then be considered as a virtual 
emitter for the outer space. 


3 Reference 2, p. 479. 
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Fic. 1. Representation of element and space potentials 
in double grid tube (longitudinal half-section). A, B, C, 
D and E, indicate successive stages in the establishment 
of multiple space charge. 


The situation may be discussed most easily by 
considering the results of experiments with a 
double grid tube (Type 24). The emitter in this 
tube is of the indirectly heated type which 
furnishes an equipotential cathode, thus avoiding 
the difficulty presented by the drop in potential 
along a filament. The 24 also has relatively large 
spacing between the second grid and plate, a 
situation which is favorable for the establishment 
of space charge in this region of the tube. 


Ill. MuttipLe SPACE CHARGE IN THE DOUBLE 
GRID TUBE (TYPE 24)—RESULTs 
AND INTERPRETATION 


Fig. 1 shows a longitudinal half-section view 
of the elements in the 24 tube. K refers to the 
cathode, G,; to the first grid, G; to the second 
and P to the plate. A set of element potentials 
under which external space charge effects are 
observed is indicated by the position of the 
points E,, E, and E, taking the cathode po- 
tential as zero. If we neglect penetration and end 
effects, we may represent qualitatively the form 
of the inter-electrode fields by drawing lines con- 
necting these points. The cylindrical arrangement 
of the elements leads to a space potential dis- 
tribution of the type shown by lines A.‘ The field 
at any point can be obtained from the slope of 
the potential curve. 

When emission occurs, all space potentials are 
depressed, the largest effect occurring in the 


* See Chaffee, Theory of Thermionic Vacuum Tubes, p. 75. 
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Fic. 2. Multiple space charge characteristic in double 
grid tube (Type 24). Plate current is plotted against 
emission (logarithmic scale) with constant grid and plate 
potentials. The three curves represent three outer grid 
potentials. 


region of weakest field (i.e., near the plate). 
With increase in emission, a cylindrical surface 
of minimum potential is established near the 
plate which travels inward toward Gs as the 
electron density grows. It is assumed that suc- 
cessive stages in this process might be indicated 
by such curves as B, C and D. 

Fig. 2 shows the variation in plate current 
with emission for several values of Ee. For con- 
venience one of these curves (E,.=4.5 volts) has 
been divided into sections. The letters which 
mark the divisions are placed to represent the 
situations indicated by the same letters in Fig. 1. 
It is characteristic of the curves shown that 
the plate current first rises to a rather sharp 
maximum level determined by the grid po- 
tentials; then, with no change in condition 
except a further increase in emission, falls, 
passes through a broad minimum and eventually 
rises slowly. While interest attaches primarily to 
the current which reaches the plate, our under- 
standing of the situations which affect the 
electron stream is materially aided by a knowl- 
edge of the grid currents. Fig. 3 gives the relative 
magnitudes of these currents expressed as frac- 
tions of the total emission. The case when 
E,=4.5 volts is taken as an example. 

The currents to all three electrodes are limited 
originally by the number of electrons present 
(i.e., by emission). Let us now follow the changes 
which occur as the emission increases. At B, 
Fig. 3, we note an increase in the fraction of the 
stream reaching Ge, and a corresponding de- 
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Fic. 3. Relative number of electrons reaching the elements 
as emission changes. 


crease in J,/I,. This is interpreted as an indica- 
tion that space charge conditions have been 
established in the outer space. If this is correct, 
limitation from this cause will retard the increase 
of current to the plate. This stage is evident in 
Fig. 2 in the decrease in slope of the plate current 
curve. 

A simple space charge effect might be expected 
to produce a saturation current level, the height 
of which would depend largely on plate potential. 
In place of this we find the current reaching a 
maximum and then declining with increased 
emission. Can this phenomenon be explained on 
the basis of space charge? 

Let us consider the influence of the surface of 
minimum potential which appears to be estab- 
lished in space 3. Only those electrons having 
sufficient velocity to move against the retarding 
field as far as this surface will contribute to the 
plate current. Those of lower speed will return 
in the direction of G2. They may strike this grid, 
or pass through its mesh toward G;. Their be- 
havior will depend on the fields encountered in 
space 2. In any case their presence near G2 will 
further lower the space potential in the vicinity. 
(Unless all returning electrons strike the grid on 
the first passage, this influence may extend to 
the potential of space 2. The effect is intensified 
if several passages through the mesh occur 
before capture by G:.) 

The maximum current to Ge is not reached 
until the emission is about 2000 microamperes. 


At 800 microamperes the curve showing J,/J, 
makes a sharp bend upward, because, it would 
seem, of the return of space charge limited 
electrons from a minimum potential region in 
space 2. Since increased emission, at least in the 
earlier part of this stage, will mean greater 
electron density in space 3, the radius of the 
outer surface of minimum potential will become 
smaller. The result is seen to be an increase first 
of Jz and then of J; at the expense of the plate 
current. 

There remains to be explained the minimum 
level of J, followed by its rise as emission be- 
comes quite large. The potential distribution in 
the tube when the plate current reaches the 
minimum stage may have the general character 
of the lines marked D in Fig. 1. Neglecting, for 
the moment the possible contribution of change 
in initial emission velocities and non-radial com- 
ponents, we may apply the three-halves power 
space charge approximation to the current from 
the minimum potential surface in space 3. Thus, 
I,=kV'/x?, where V is the potential of the plate 
with reference to the minimum, and x the 
distance from the plate at which the minimum 
is established by space charge. 

Let us now consider the influence on this 
current of the potential minimum produced in 
space 2 by the increase in density of the emitted 
stream. Electrons which now pass through G, 
will enter space 3 with higher average velocities. 
It should be recalled that the current density at 
any point is given by i=pv where p is the 
electron density. If the number entering space 3 
remains constant (due to limitation in space 2) 
this can result only in a decrease in space charge 
near Ge. Hence the potential minimum region 
should move outward toward the plate, and with 
decreasing x an increase in plate current should 
follow. The writer’s interpretation of this change 
is indicated at E in Fig. 1. The actual effect on 
the plate current is shown at E in Fig. 2. 

The influence of the potential of the second 
grid is clearly shown in Fig. 2. The curves for 
various potentials branch off from the early part 
of the curve which is, within experimental error, 
common. It seems significant that with in- 
creasing E» the peak shows only a slight shift 
toward higher total emission, while the amplitude 
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of the effect, i.e., the spread between maximum 
and minimum values of current, grows larger. 
Both of these trends are in line with the explana- 
tion given above. For example, a higher po- 
tential difference between Gz and P will require a 
greater stream density to produce space charge 
limitation. Once established, however, a larger 
fraction of the total will finally reach Gz. This 
has been verified by the comparison of the curves 
in Fig. 3 with similar curves obtained with 
different values of Ee. 


IV. Errect or CATHODE SPACE CHARGE 


It has been assumed in the foregoing that 
space charge limitation in the space between the 
emitter and first grid was not present when the 
phenomenon was observed. This has been shown 
experimentally to be true in the range thus far 
considered. It is reasonable to inquire into the 
influence of E; on the observed behavior of the 
electron stream, and in particular to test the 
effect of inner space charge. 

When £; is varied, two significant effects are 
observed (see Fig. 4): (1) For lower E, the peak 
becomes larger and suffers a displacement in 
the direction of the lower total emission ; (2) the 
final rise of plate current occurs earlier (at 
lower emission values) for low than for high E. 
Any change which results in more readily estab- 
lishing space charge in space 2 will cause the 
second of the trends. This may also be said for 
space charge surrounding the emitter. This 
probably accounts for the large difference ob- 
served between the 6 volt and the 16 volt 
curves. 
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Fic. 4. Multiple space charge characteristic in double grid 
tube. Effect of Fy. 
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The displacement of the peak presents a more 
difficult problem. The explanation which seems 
probable at the present is that non-radial com- 
ponents of velocity become more prominent for 
high grid potentials. An increase in the magni- 
tude of such components would increase space 
charge without contributing to the plate current. 


V. DIscussION * 


Two aspects of the multiple space charge 
characteristics have particularly interested the 
writer and it seems pertinent to call attention to 
them here. It has been shown theoretically and 
supported experimentally that the fluctuation 
level (E?) of the normal shot effect is reduced 
under space charge by a factor (@/,/0I.)? where 
I, is the emitted current and J, that which 
reaches the collector.5»® For two points on the 
characteristic reported here this factor is zero. The 
same consideration may be extended to cover 
small changes in emission due to variations in 
the heating current or changes in the emitting 
surface.’ As long as these are of the relatively 
small amplitude usually encountered the effect in 
the plate circuit is very small indeed. 

It has further been noted that the electron 
stream reaching the outer space is pure (i.e., free 
from positive ions). One is therefore able to 
establish space charge conditions in this region 
without introducing current fluctuations of the ab- 
normal type described by Kozanowski, Williams 
and the writer.® ° 

The fluctuation level in the output circuit of 
a tube operated in this manner may thus be 
reduced practically to that resulting from thermal 
agitation. If, therefore, a control grid be inserted 
between the virtual emitter and plate, the tube 
may be operated as an effective triode with the 
unique property that the plate current is inde- 
pendent both of statistical and ‘‘functional” 
Variations in emission. 


°F. B. Llewellyn, Proc. I. R. E. 18, 243 (1930). 

6 See reference 2, p. 115. 

7 E.g., the “‘flicker’’ effect which becomes prominent at 
low frequencies. Johnson, Phys. Rev. 26, 71 (1925); 
Schottky, Phys. Rev. 28, 74 (1926). Pearson, Physics 5, 
233 (1934) has demonstrated the effectiveness of space 
charge in eliminating these fluctuctions in the space 
current. 

8 Kozanowski and Williams, Phys. Rev. 36, 1314 (1930). 

® Thatcher and Williams, Phys. Rev. 39, 479-484 (1932). 





MULTIPLE SPACE CHARGE 


It is well known that one of the limits on the 
useful gain of vacuum tube amplifiers is the mag- 
nitude of these fluctuations in early stages— 
particularly the first. The operation of the first 
the manner described should make 
possible an increased signal to noise ratio in the 
input to the second tube. 

The properties of multiple space charge seem 
of particular interest in d.c. amplifier applica- 
tions of the electrometer type. In this connection 


tube in 


several schemes have been devised to minimize 
the effect of filament variation.'® As far as the 
writer is aware, however, all these have involved 
the matching of tubes or setting up bridge 
circuits by which variations in one branch are 
compensated by those in another. The use of the 
multiple space charge characteristic should re- 
duce the complexity of these operations and 
increase sensitivity by elimination of one of the 
main causes of instability. 


VI. RESULTS WITH TRIPLE GRID TUBE 
(TYPE 78) 


A preliminary test of this method of operation 
has been made in a study of the phenomenon in 
a triple grid tube (Type 78). The tube chosen is 
one with external connections to each grid. 
Field conditions similar to those in the preceding 
case were established surrounding the emitter. 
The third grid was used as a control element. 
Fig. 5 shows the plate current emission charac- 
teristics which were obtained with £; = +9 volts, 
E,=+4.5 volts and E,=+9 volts. 

Attention is directed to the displacement of 
the initial peak as E; is varied. This peak is 
reached earlier as E; is made more negative, 
with consequent increase in outer space charge 
and reduction of plate current. The relatively 

’ See, for example: Soller, Rev. Sci. Inst. 3, 416 (1932); 
Du Bridge and Brown, Rev. Sc:. Inst. 4, 532 (1933); 
peda Van Voorhis, Rev. Sci. Inst. 5, 244 (1934); 
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minimum occurs at about the same 
position on each curve. With 4500 microamperes 
emission, the equilibrium potential of Gs; was 
determined. This is the potential which the grid 
will assume when left free from connection with 
other parts of the system. If a bias of this amount 
is applied, no change in current will be ob- 
served. 


broad 


For the tube considered here this value was 
found to be —0.21 volt, at which potential a 
stream of 23.8 microamperes reaches the plate. 
The transconductance in this region was found 
to be about 200 micromhos. With tubes designed 
to provide more effective control on G; this figure 
should be materially exceeded. Improvement can 
particularly be expected in the increase in input 
resistance which will result from bringing the 
lead from this grid to the cap, rather than 
through the pinch seal to the base of the tube. 
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Some Electrical Properties of Ceresin Wax 


A. Wixstrom, Cornell University, Ithaca, N. Y. 
(Received October 26, 1934) 


A lower limit for the resistivity of commercial ceresin is 
set at 1.76X10'* ohm-cm if the voltage gradient is 2000 
volts per cm the thickness of the sample is 0.214 cm and 
the temperature is 23°C. There is a comparatively large 
scattering of time functions used in the equation for the 
discharge current if expressed in accordance with Wagner’s 
theory. The charge Q “‘stored” in ceresin is a function of 
time. At any time after the application of the voltage, this 


charge is proportional to the voltage. PE was measured 
as a function of t. The increase is found to be more rapid 
than is claimed by other investigators. No variation of 
P/E with E was found. The results, however, are not quite 
conclusive, because of the inherent weakness of the method 
used in determining the polarization. The resistance R as 
defined by Richardson is calculated and found not to be 
constant for ceresin but to vary both with time and voltage. 





STUDY has here been undertaken of 

ceresin wax when subjected to a constant 
voltage, in order to determine some of its 
electrical properties. 

The sample used in this investigation was 
obtained from a block of commercial ceresin, no 
attempt being made to purify the product. The 
reason was that this wax, which has a yellowish 
color, is used today to a very great extent in 
research laboratories as an insulator. Very little 
is, however, known with respect to its behavior 
when subjected to a constant voltage. Troubles 
have often been encountered due to its soakage 
current which prevents the conduction current 
from getting reasonably constant over a period 
of several hours after the voltage has been 
applied. 


APPARATUS AND EXPERIMENTAL PROCEDURE 


For this investigation a test-condenser was 
made in the following manner. A circular brass 
plate (diameter=3 inches, thickness=}3 inch) 
was set into a lathe and the surface was carefully 
turned. The chuck was then unscrewed and laid 
flat on a table with the brass plate up. A paper 
form was made to fit the circumference of the 
plate and after the plate had been heated with 
the flame from a Bunsen burner to approximately 
80°C, melted ceresin was poured onto the brass 
plate. The paper form prevented the wax from 
flowing off the plate. After 24 hours the chuck 
was screwed back on the lathe and the ceresin 
was turned down to a thickness of 0.0842 inch 
(2.14 mm). The surface of the wax was then 
covered with several layers of gold leaf such as 
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bookbinders use. While the condenser was still 
in the lathe a circular cut was made through the 
gold to form a guard ring. 

A simplified wiring diagram is shown in Fig. 1. 
E is a Compton electrometer of low capacity, 
adjusted by means of negative control to a 
sensitivity of approximately 10,000 mm per volt, 
which could be checked with the potentiometer 
P,. The test condenser C;, shown with grounded 
guard ring, has one plate connected with the 
electrometer and the other plate connected with 
an electromagnetic switch S,, which made it 
possible to connect C; either with the ground or 
with one terminal of the high tension battery B, 
the other terminal of which is grounded. C, is a 
cylindrical condenser (see Fig. 2), the capacity 
of which is calculated from its dimensions using 
the formula: C=L/21n (d;,/dz). The quantities 
L, d, and dz were chosen so as to make C= 25 cm. 





Fic. 1. Diagram of connections for the determination of 
the total current 7, or discharge current iq. 





Zz 











Fic. 2. Cylindrical condenser with guard rings. 
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Fic. 3. Diagram of connections for the determination of 
the polarization P. 


The dotted lines in Figs. 1 and 3 show to what 
extent the apparatus was shielded. Besides this 
shielding the switches S, and S2 were in a separate 
grounded metal box, as these switches were 
operated electromagnetically. The shielding box 
containing the specimen had for drying purpose 
a large surface dish containing calcium chloride. 

The time recording was done by means of a 
chronograph manufactured by. Wm. Gaertner 
and Company, Chicago. Before being put into 
service this instrument was adjusted and checked 
by means of a synchronous clock motor. The 
maximum error was +0.005 of a second per 
minute. 

In this investigation the following quantities 
were measured : 

(1) The total current 7, vs. time ¢ for different 
values of voltage E of the battery B (see Fig. 1). 

(2) The discharge current ig vs. time ¢ for 
different values of the time 7 during which the 
voltage E was applied to the dielectric. 

(3) The polarization P vs. charging time T. 
One such curve obtained for each voltage E. 


(1) The experimental determination of i; 


Before any measurements of the total current 
i; were made, the condenser C; was left short- 
circuited until completely discharged. This was 
done by throwing the switch S; into position 
No. 1 and closing K, for several days. The 
condenser was considered completely discharged, 
when the insulated electrometer quadrants upon 
opening K, would not take on a charge at a 
greater rate than 5X10~-'® coulomb per second. 
When this was the case the total current i; was 
determined as a function of time in the following 
manner. * 

At a certain instant of time ‘=0 the switch 
S; was moved into the position No. 2 (see Fig. 1), 


while K, and Kz were closed. The current, which 
will then flow, consists of a transient current 
superimposed upon the total current 7,. Because 
of the magnitude of the constants of this circuit 
the transient current subsided within a fraction 
of the first second. No attempt was made to 
measure it. 

Whenever 7, was to be determined K, was 
first opened. The electrometer then drifted and, 
when passing a certain division N upon the scale, 
the time ¢=¢, was recorded. The electrometer 
was now brought back to zero and kept there by 
adjusting the potentiometer P:. Towards the 
end of this determination the continuous adjust- 
ment of P2 was stopped, resulting in a drift of 
the electrometer. When again the Nth division 
on the scale was passed, the time t=f was 
recorded. The switch K, was now closed and as 
soon as the electrometer was back to zero, 
another determination of i; could be made. 

The voltage change AE of the potentiometer 
P: was measured with a Leeds and Northrup 
potentiometer. If the capacity of the standard 
condenser was C,, the average value of 7, between 
t=t, and t=#, can be calculated from the equation 


i, =AQ/At=C,AE/(tp—h). 


This average value was associated with the 
instant of time ¢=(4+42)/2. It is fully realized 
that this is not strictly correct. The function 
“current vs. time’’ may even for small values of 
At, deviate appreciably from a straight line. A 
method of measuring instantaneous values of 
rapidly decreasing currents is described in a 
separate paper by the author.’ For values of 
t<1 minute this method was found desirable 
and generally used. 


(2) The experimental determination of i, 


After a certain length of time t=T7, the test 
condenser C; was completely discharged by 
closing K, and throwing the switch S, into 
position No. 1 for a sufficient length of time. 
Immediately after such a condenser is short- 
circuited a current will flow which consists of a 
transient current superimposed upon the dis- 
charge current i,. Again, no attempt was made 
to measure the transient. The discharge current 


1A. Wikstrom, Rev. Sci. Inst. 4, 612 (1933). 
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i, Was measured in the same way as 7% and 
readings were continually taken until the values 
were of the order of 5 10~-'* ampere or less. 

It was mentioned above that the condenser 
was considered to be completely discharged 
when, upon opening K,, the insulated pair of 
quadrants did not charge up at a greater rate 
than 510-'® coulomb per second. Even when 
the condenser was left short-circuited for more 
than nine days this rate of charge varied in an 
irregular manner from day to day between 
1x10-'® and 5X10-'* coulomb per second. Re- 
versing the needle potential of the electrometer 
had no effect upon either the magnitude or the 
sign of the rate of charge accumulated which 
was always found to be of positive sign. Thus it 
was not due to a leak of charge between the 
needle and the insulated pair of quadrants. 
When the electrometer was disconnected from 
the rest of the circuit, it took on charge at a rate 
of less than 10-'? coulomb per second. This 
indicates that the charge was not released from 
the electrometer insulators. It was only when C, 
was connected with the electrometer that the 
rate increased to a value within the 
mentioned limits. 

With the same polarity of the needle the 
electrometer, when insulated, always drifted in 
the same direction. The readings obtained were 
therefore either too large or too small, depending 
upon the direction of the current measured. 
The values of 7, were always so large as not to 
be affected much by this unidirectional error. 
To eliminate it as far as possible it was, therefore, 
necessary to repeat every 7, curve under identical 
conditions, with the exception of E being re- 
_ versed. The correct value of 7, was then taken to 
be the arithmetic mean of corresponding values 
of the two current curves obtained in these two 
sets of readings. 


above 


(3) The experimental determination of P 


Richardson's? definition of polarization was 
used and except for the mechanical manipulation 
the method of measuring it was identical with 
his method. A simplified diagram is shown in 
Fig. 3. Before any reading was taken the con- 
denser C, was completely discharged as in the 


2S. W. Richardson, Proc. Roy. Soc. A92, 101 (1915-16). 
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case of the current measurements. At an arbi- 
trary time ¢=0 5S, was thrown into position 
No. 2, K, of course being closed. If now K; was 
opened, the electrometer E began to drift. The 
direction of this drift was noticed, whereupon 
K, was again closed. We shall call this direction 
of drift the positive direction. At a certain 
instant of time t=7, at which the polarization 
was to be determined, S2. was changed into 
position No. 2, thus applying to C; a voltage 
E’ <E. Immediately thereafter K,; was opened 
and the direction of the drift of the electrometer 
was again observed. If the polarization voltage 
is denoted by P, then a positive drift indicates 
that E’>P, whereas a negative drift shows that 
E’ <P. This procedure was repeated, beginning 
every time with a completely discharged con- 
denser, until E’=P was obtained, as was indi- 
cated by no drift of the electrometer for an 
appreciable time after that K; was opened. 

The switches S,; and S. were operated by 
electromagnets. A double pole knife switch closed 
the circuit containing the electromagnet which 
operated S; and at the same time recorded 
through a separate circuit the time by means of 
the chronograph. 

The grounding key K, was operated by means 
of a lever and a string. The string, which when 
pulled opened K,, was fastened to a single pole 
knife switch of special design. When operated, 
this switch first, through an electric circuit, 
threw 5S: into position No. 2, and then K, was 
immediately opened by means of the string. It 
was mechanically impossible for K, to open until 
Se made contact in position No. 2. 

The battery B’ consisted of 1.5 volt dry cells, 
thus making a fine adjustment of E’ possible. 
A potentiometer, not shown in Fig. 3, was also 
included for the same purpose. 

In order to minimize disturbances due to 
temperature variations, the research room was 
kept at 23°C with a maximum variation of 
+0.5°C. 


EXPERIMENTAL RESULTS 


A preliminary test showed that the volume 
resistivity of paraffin is much higher*than that 
of ceresin. As the sensitivity of the electrometer 
used in this work was not high enough to make 





PROPERTIES 


any comparative measurements, the relative 
magnitude of the resistivities of ceresin. and 
paraffin cannot be given. Strictly speaking, there 
is no one definition of “‘resistivity’’ which can be 
used for dielectrics in general. Various attempts 
have been made to define some such constant 
that will make it possible to apply Ohm’s law 
to dielectrics. To the author’s knowledge, so far 
no definition has proved successful. A typical 
example is dealt with later in this paper. 


Data directly measured 


The results of direct measurements are shown 
in Figs. 4, 5 and 6. These curves are taken as 
representative of a comparatively large amount 
of data accumulated during the course of the 
investigation. Plots of 7, against the charging 
time are shown in Fig. 4 on log-log paper where 
each curve corresponds to a different constant 
potential applied to the condenser. These curves 
are fairly straight lines from ¢=0.1 minute to 
t=1 minute. 

In general the total current was only measured 
for 25 minutes and did not become reasonably 
constant during this time. In one case the current 
measurements were continued for 1360 minutes. 


Fic. 4. The total current i, vs. time. 
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Fic. 5. The discharge current i, vs. time for E=651 volts 
for four different charging times 7. 


The following calculations show how the re- 
sistivity p was determined for this particular 
voltage : 

p=EXxD?/i,X4d, 

E=425 volts, 

t,=4.91X10-" amp., 

D=7.436 cm, 

d=0.214 cm. 


Hence p=1.76X10'® ohm-cm. 


This calculation is made in order to set a 
lower limit to the resistivity of ceresin wax under 
the following conditions: A voltage gradient of 
nearly 2000 volts per cm; a thickness of 0.214 
cm and a temperature of 23°C. The value for p 
found by Guében?® is 3 x 10'8 ohm-cm. 

The discharge currents 7, were determined for 
four different values of charging time 7°, namely: 


3G. Guében, Zeits. f. Physik 74, 714 (1932). 





A. WIKSTROM 


Pexpressed in percent of E vs time. 







































































70 2 4 i 2 22  & 
t minutes 


ft ¢€ § 


Fic. 6. The polarization P (expressed in percent of E) 
vs. time. 


1, 5, 10 and 25 minutes. The curves of discharge 
current vs. time for E=651 volts are shown in 
Fig. 5 plotted on log-log paper and are seen to be 
slightly bent lines. As 74 decreases faster than a 
straight line it is possible that for very high 
values of ¢, the function i(t)g will gradually 
become an exponential function. It was not 
possible to determine whether or not this is true, 
as it occurs, if at all, at such small values of i, 
that the accuracy of the measurements would 
not suffice. 

The polarization P expressed in percent of E 
is given as a function of charging time in Fig. 6. 
From the author’s observations there is no 
evidence of a variation of P/E with the applied 
voltage. 


Data calculated from measurements 


From curves shown in Fig. 5 and many similar 
curves obtained at other voltages, the residual 
charge Q in the dielectric was determined. The 
discharge current for a given time 7 of charging, 
is a function of the time during which the dis- 
charge takes place and if this current is expressed 
by the symbol i(f)a, the charge Qr, which is 























Fic. 7. Qr vs. t. The charge Qr “‘stored”’ in the dielectric 
vs. time. Each curve is for a constant potential E applied 
to the test condenser. 
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Fic. 8. The charging current dQ/dt vs. time. 


accumulated in the dielectric during the charging 
time 7, can be written in the following manner: 


Or= f i(t)adt. 


As it is impractical to evaluate Q7 in a theo- 
retical way, the discharge current i(¢)4 was 
plotted to large scale on graph paper, and the 
area under the curve was determined with a 
planimeter. The values of Q7 determined in this 
manner are plotted in Fig. 7 against T as abscissa. 
For each curve the voltage is constant. 

Plotting Q7 against E for constant values of T 
gives the linear relations shown in Fig. 9. 
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Fic. 9. The charge Q7 ‘“‘stored” in the dielectric vs. the 
applied potential E. Each curve is for a constant charging 
time T. 


¢ minutes 


Fic. 10. R vs. t. The resistance R calculated from 
Richardson’s equation as a function of time. Each curve 
is for a constant potential E. 


The derivative of Qr with respect to the 
charging time may be taken to be the charging 
current often spoken of as the soakage current. 
From the curves shown in Fig. 7 dQ/dt is plotted 
against charging time as abscissa. Thus a family 
of curves is obtained as is presented in Fig. 8. 
For each such curve the voltage is constant. 

Since the author used Richardson’s definition 
of polarization it was of some interest to compute 


the resistance R as defined by Richardson‘ in 
his equation for dielectrics 
1,=dQ/dt+(E—P)/R. 


Fig. 10 shows this resistance vs. time for different 
voltages. 


Discussion of results 


The discharge curves shown in Fig. 5 are 
slightly bent lines. This, if compared with 
Wagner’s® theory, calls for a small value of 6 in 
the equation 


etc 
I(t) = (ekobEo/4200v 7) f ee 22— 2z—te— */00d 2, 


which means that there is a comparatively large 
scattering of the time functions, or that the time 
constants do not group themselves very closely 
round the most probable value In 4p. 

In Fig. 6 is shown the variation of P/E with 
time. Owing to the rapid increase in the polariza- 
tion with time and the inherent inaccuracy of 
the method used in determining P, we should 
not give too much weight to the statement that 
P/E is independent of E. The results shown in 
Fig. 6 are averages of many readings. The 
previous condition is, in these measurements, of 
the greatest importance and care was therefore 
exercised so as to leave the dielectric short- 
circuited for a long time before a new reading 
was taken. 

While Courtener and Ignatowitsch* found that 
for ebonite P/E is independent of E, Tutschke- 
witsch’ concludes that this is not true for paraffin. 

A direct comparison of data for ceresin has 
been made possible by the work of Nasledow 
and Scharawsky.* Comparing the polarization 
found by these investigators with the author’s 
values at corresponding values of time it is found 
that the author’s values are higher. This differ- 
ence may be due to different composition, 
temperature and thickness of the samples. Thus 
the difference in temperature of almost 6°C 


4S. W. Richardson, Proc. Roy. Soc. A92, 41 (1915-16). 

5K. W. Wagner, Ann. d. Physik (4) 40, 817 (1913). 

6 A. Courtener and J. Ignatowitsch, Ann. d. Physik [5] 
10, 681 (1931). 

om M. Tutschkewitsch, Ann. d. Physik [5] 6, 622 
(1930). 

8D. Nasledow and P. Scharawsky, Ann. d. Physik [5] 
5, 429 (1930). 
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would, on the basis of the dipole theory as 
modified by Whitehead,’ account for the quicker 
polarization obtained by the author. An increase 
in temperature decreases the coefficient of vis- 
cosity and the dipoles therefore line up faster. 

The curves shown in Fig. 7 are similar to 
values obtained by Courtener and Ignatowitsch® 
and Nasledow and Scharawsky’® all indicating 
that the charge accumulated in a dielectric 
approaches a maximum value of charge that 
can be stored in the dielectric itself. For higher 
values of voltage gradient it takes the same time 
to reach a value which is within a certain percent 
of Quan. 

The investigators referred to above have all 
studied samples of different thickness. It is 
conceivable that the thickness may have various 
effects upon the characteristics of a dielectric. 
Thus the rapidity by which the polarization 
takes place may possibly change contrary to the 
theory of similarity suggested by B. B. Kurt- 
schatow and J. Kurtschatow.!® 


S. Whitehead, Phil. Mag. [7] 9, 865 (1930). 
B. B. Kurtschatow and J. Kurtschatow, J. d. 
phyvsik. chem. Gesell. 61. 321 (1929). 
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SADRON 


The purpose of this investigation was to study 
ceresin such as it is believed to be found in use 
for insulation in research laboratories of today. 
The characteristics can probably be changed to 
a very great extent, depending upon the way 
the wax is treated while the condenser is made. 
In order to permit results to be compared 
successfully, a standard treatment should be 
adhered to, such treatment including the heating 
of the wax to the melting point for a definite 
time in an evacuated chamber and possibly in a 
strong electrostatic field. It is believed that the 
electric field will partly free the wax from 
heterogeneous parts which are usually present 
in the form of ions from impurities. 

In conclusion the author wishes to thank Dr. 
L. L. Barnes and Professor E. H. Kennard for 
many helpful suggestions and much valuable 
criticism in connection with this investigation. 
He is further much indebted to the members of 
the Department of Physics, and in particular to 
Professor Ernest Merritt and Professor F. K. 
Richtmyer, who have permitted him to make use 
of the equipment and research room which were 
needed to carry out this work. 
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An Optical Method for Measuring the Distribution of Velocity Gradients in a Two- 
Dimensional Flow 


E. D. Atcock anp C. L. Sapron, California Institute of Technology, Pasadena, California 


(Received November 15, 1934) 


I. ANISOTROPY OF 


AXWELL that some _ very 
viscous fluids (Canada balsam) when in 
motion present the phenomenon of double re- 
fraction. Numerous scientists have studied the 
phenomenon from experimental as well as theo- 
retical points of view, the first real work being 
done by Kundt in 1881. The result of these 
investigations is that all colloidal suspensions and 
a large number of pure liquids show a consider- 
able double refraction when in motion. 
In the case of colloidal suspensions, theories 
have been advanced by several authors, not 


Liguips IN LAMINAR FLOW 


discovered 


always with great success, since the experimental 
work shows that the effect is rather complicated. 
For pure liquids we have only a theory of Raman 
and Krishnan which is still open to discussion. 
Although the theory in this case is more difficult 
than in the case of suspensions, the experimental 
data indicate a simpler effect. When we are 
Cealing with pure liquids we treat them in the 
following manner: 

Let us consider a parallel or laminar two- 
dimensional flow along a solid boundary and 
suppose there is a constant gradient of velocity 
perpendicular to the boundary. Such a flow will 
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Fic. 1. Experimental arrangement. S, arc; L,L.L;, lenses; N,N2, nicol prisms; C;, gypsum wave 
plate; F, spectrograph slit; P, prism; P’, photographic plate. 


have shearing stresses parallel to the boundary. 
These stresses may be resolved into a tension and 
a compression perpendicular to each other and 
making an angle of 45° with the boundary. The 
effect of these stresses on the optical properties 
of the liquid is analogous to the effect produced 
in photoelasticity when a solid material (Bake- 
lite, glass, etc.) is subjected to stresses. If we 
now send a beam of polarized light. through the 
liquid perpendicular to the plane of the flow, the 
electromagnetic vector is resolved into two com- 
ponent vectors along the axes of tension and 
compression. The two components travel with 
different velocities due to the double refracting 
properties of the fluid. If the light path is of unit 
length, the separation of the two vibrations is 
given by 
D=n,—Ms, 


where 2, and m2 are the indices of refraction in 
the two principle directions. The experiment 
(and the theory) show that the difference of the 
indices (%,;—2) is proportional to the stresses 
and, since a parallel flow was assumed, we have 
Ou 
D=yK—, 


oy 
where K is the so-called Maxwell constant which 
depends upon the molecular properties of the 
fluid and yu is the coefficient of viscosity. 

If the velocity gradient is a function of y ina 
given cross section, D will be proportional to 
the same function of y. If our flow is curved, the 
double refraction at any point is given by 


Ou 
D= uk —u RI 
oy 


where R is the radius of curvature of the flow 
at the given point. 


Il. APPLICATIONS 


In all cases where the radius of curvature of 
the flow is large so that the term uw R may be 
neglected, the above phenomenon furnishes us 
with a simple method of measuring the distribu- 
tion of velocity gradients or of velocity without 
disturbing the flow. The problem is reduced to 
measuring the double refraction D and Maxwell's 
constant K. From these two determinations we 
obtain directly the velocity gradient and by a 
simple integration, the velocity. It is a simple 
problem in crystalline optics to determine the 
double refraction. We used a method already 
used by several other experimenters. The double 
refraction produced by the liquid is usually very 
small, therefore we use it to modify the double 
refraction of a gypsum wave plate. By measuring 
the shift in the interference fringe of the wave 
plate with a spectrometer we obtain the double 
refraction of the flowing liquid. 


First application: Parallel channel 


In order to check the value of our method we 
applied it to the well-known case of flow between 
two infinite planes. We used a channel of 
rectangular cross section with a very large aspect 
ratio. (See Fig. 1.) In this channel we circulated 
sesame oil with an average velocity of about 1 
m ‘sec. by means of a pump loaned to us through 
the courtesy of the Byron Jackson Company. 
We measured this average velocity by taking 
motion pictures of the falling oil in the gauge 
and of a stopwatch simultaneously. Theoreti- 
cally the flow in a cross section of this type gives 
us a linear distribution of velocity gradient from 


one wall to the other. Fig. 2 shows our experi- 
mental results, which are in quite good agree- 
ment with the theory. Because of strains in the 
glass side walls of our channel the fringes are 
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Fic. 2. Distribution of du/dy as a function of y in the section Oy. 


Fic. 3. Distribution of du/dy in cross section of diverging channel. 


Fic. 4. Distribution of u (velocity) in cross section of diverging channel. 





DISTRIBUTION 


oO 


300 400 Se 
7) 100 200 > ale 


Fic. 5. Double refraction in fluid between 
rotating cylinders. 


not straight ; however, this cancels out since we 
only measure the shift of the fringe. By com- 
bining the theory and experiment we were able 
té calculate a value of Maxwell’s constant of 
9x10-" c.g.s. Using an aspect ratio of 1/30, we 
computed the end effects. This calculation 
showed that we were introducing an error of not 
over 3 percent by neglecting them. The esti- 
mated accuracy of the method is about eight 
percent. 


Second application: Diverging channel 


We then applied the method to the problem 
of the flow in a diverging channel. The results 
are given in Figs. 3 and 4. Since the quantity of 
fluid flowing per unit time through any section 
is constant, the integral of the velocity across 
any section should also be constant. We have 
indicated in Fig. 4’ the results of such an inte- 


gration for several sections and they check very 
well. 


Third application: Rotating cylinders 
For the purpose of checking our value of 
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Maxwell's constant obtained in our first experi- 
ment and to investigate what happens when the 
flow becomes turbulent, we took the case of flow 
between two concentric cylinders. The inner 
cylinder rotates while the outer remains station- 
ary. We had several cylinders made which 
enabled us to measure the distributions in five 
different gaps. The results are shown in Fig. 5. 
The average value of Maxwell’s constant calcu- 
lated from these results is 8.5 X10-" c.g.s. which 
checks very well with that obtained in our first 
experiment. Using sesame oil we were not able 
to attain a high enough velocity to have turbu- 
lent flow. In order to have a control on the 
temperature of the oil we installed a water jacket 
on the outer cylinder. 


CONCLUSIONS 


We think that our method is applicable to 
many problems, it being subject only to the 
limitation that the flow be two dimensional. 
We see at present two main fields of application. 

(1) Scientific investigation of the boundary 
layer. If sufficient care is taken we estimate that 
it is possible to measure the value of the velocity 
gradient to within 0.1 mm of the boundary. 
Therefore we can obtain the velocity distribution 
in this layer without disturbing the flow with 
any apparatus. Also since we can apply the 
method to curved flows where the term u/R can 
be neglected we can investigate several problems 
connected with separation. 

(2) Experimental investigation of turbulent 
distributions. We cannot give any details at 
present on this problem, our work not being far 
enough advanced. 

Finally we wish to point out that this method 
could be easily applied to the problem of the 
flow around a rotating system. Instead of using 
a continuous source in this case, we would have 
to use a stroboscopic device. 
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Magnetic Properties and Orientation of Ferromagnetic Particles' 


C. W. Davis, U. S. Bureau of Mines, Pittsburgh, Pennsylvania 
(Received October 19, 1934) 


Results of preliminary tests by direct methods show definitely that alignment of particles in 
a magnetic field increases their magnetization up to saturation and augments their remanence, 
but that coercive force is not noticeably affected. 


INTRODUCTION 


HILE engaged in determining magnetic 

properties of mineral powders, incidental 
to the study of mineral physics being conducted 
by the metallurgical division of the U. S. Bureau 
of Mines,* the writer observed that magnetic 
particles of a definite size and packing density, 
when oriented by placing in a strong magnetic 
field and then demagnetized, gave higher values 
for remanence and induction at moderate field 
(45 than did undirected 
particles under similar conditions. After a satis- 
factory modified isthmus method had been de- 
veloped‘ for determining magnetic properties of 
powders at field strengths up to about 3600 
oersteds, the effect of particle orientation on 


strengths oersteds) 


these properties was studied in more detail. 


EXPERIMENTAL 


Natural and artificial magnetic oxides of iron, 
the purification and preparation of which are 
described elsewhere,’ were studied. 

Natural Fe;O,, a sample of well-crystallized 
magnetite from Mineville, New York, was sized 
between 60 and 250 mesh and found to contain 
96.97 percent Fe;O, (calculated from the FeO 
content), 0.83 percent Fe2O;, 1.38 percent SiO» 
and 0.57 percent TiQs. 


! Published by permission of the Director, U. S. Bureau 
of Mines. (Not subject to copyright.) 

* Associate chemist, U. S. Bureau of Mines. 

®*R. S. Dean, and others, Progress Reports— Metal, 
lurgical Division; 10. Mineral Physics Studies, Report of 
Investigations, 3268, U. S. Bureau of Mines (in course of 
publication). 

*V.H. Gottschalk and C. W. Davis, An Apparatus for 
Determining the Magnetic Constants of Mineral Powders, 
Report of Investigations, 3268, U.S. Bureau of Mines (in 
course of publication). 

°C. W. Davis, Magnetic Properties of Mineral Powders 
and Their Significance, Report of Investigations, 3268, 
U.S. Bureau of Mines (in course of publication). 
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The artificial oxide produced by the partial 
gaseous reduction of hematite from Missouri was 
porous and impure, being a mixture of 60.07 
percent Fe;0O, (calculated from the FeO con- 
tent), 29.5 percent FesO;, 0.58 percent MgO and 
7.93 percent SiO». It was selected because of the 
exceptionally high coercive force exhibited by 
material prepared in this way® and was crushed 
and sized between 80 and 250 mesh. The minus 
250-mesh particles from each sample were dis- 
carded because it had been learned that the 
effect of size on the magnetic properties of 
powders became noticeable at smaller sizes.’ 

Magnetic measurements were made with the 
previously mentioned isthmus apparatus, using 
the point-by-point method for hysteretic de- 
terminations and taking care that all conditions 
except particle orientation remained constant 
for each pair of tests. This was accomplished by 
using a single sample of such a size that all of 
it could be packed into a sample capsule of 
fixed, known dimensions for both the random and 
regulated particle orientation; that is, the mag- 
netic properties of the same material were de- 
termined with two different particle orientations 
at the same packing density and at the same 
field strength. 

The observation of Weiss* that large mineral 
crystals had preferred direction of magnetization 
was extended te small crystals of magnetite ina 
very simple way. Controlled of 
mineral grains was accomplished by the gradual 
addition of the particles while the sample capsule, 
placed between the poles of the magnet of the 
isthmus permeameter, was agitated by tapping. 


orientation 


®V. H. Gottschalk and C. W. Davis, A Magnetic 
Material of High Coercive Force, Nature 132, 513 (1933). 

7V. H. Gottschalk and F. S. Wartman, Magnetization 
Curves for Magnetite Powders, Report of Investigations, 


3268, U. S. Bureau of Mines (in course of publication). 


* Pierre Weiss, Thesis, Paris, 1896. 
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This method may or may not orient crystals in 
any definite crystallographic arrangement, but 
it must turn particles so that the direction of 
their resultant polarity approaches parallelism 
with the magnetic field in which they become 
oriented, and by rotation of the container, in 
which mineral particles of more than one axis of 
maximum magnetization are suspended in a 
liquid gelatin solution, a very definite mineral 
orientation should result.’ 

The use of a binder to fix the position of the 
particles was found to be unnecessary for our 
tests. 

Demagnetization was effected by placing the 
sample within its capsule in a slowly reversing 
magnetic field that was gradually reduced to 
zero. 

The differences between the magnetic prop- 
erties of the samples in which the particles were 
aligned and those with random particle arrange- 
ment are shown by hysteresis loops. 


* A. Goetz, A. B. Focke and A. Faessler, Phys. Rev. 39, 

(1932); A. Goetz and A. Faessler, Phys. Rev. 
, 1053A (1932); K. S. Krishnan, Phys. Rev. 45, 115 
1934); A. Goetz, Phys. Rev. 45, 282 (1934). 
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RESULTs!” 


Fig. 1 shows that the magnetization of aligned 
particles of magnetite is considerably higher at 
moderate field strengths than that of random 
particles, the difference increasing up to about 
230 oersteds and then decreasing as the field 
increases so that each sample approaches a 47/ 
value of 2850 gauss at a field strength of about 
2700 oersteds. Both the percentage difference in 
magnetization and the difference (u;—.) in 
induction per unit field strength seem to be 
greatest for the two cases at the lowest field 
strength used, but it may well be that the 
maximum is really reached at 19 oersteds, the 
coercive force of the sample. 

The percentage difference in 47/’s and the 
difference in y's for reduced hematite samples 
increases as shown in Fig. 2 so that a maximum 
is reached at a field strength between 230 and 
254 oersteds, a value slightly higher than the 
coercive force of 230. 

The figures show that the coercive force is 
practically independent of particle orientation, 


'° Data from which the curves in Figs. 1 and 2 were 


prepared show clearly the results and relationships discussed. 
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Fic. 1. Hysteresis loops and magnetization curves for magnetite. 
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Fic. 2. Hysteresis loops and magnetization curves for reduced hematite. 


the small apparent difference that may be ob- 
tained by large-scale plots being less than the 
probable experimental error for the determi- 
nation. 

The remanence, however, is seen to be sub- 
stantially greater for the tests in which particles 
were lined up in a magnetic field, 165 as against 
137 for magnetite and 610 as against 495 for 
partly reduced hematite. 

The slight irregularities in the permeability 
values may or may not be accidental. 


DISCUSSION 


The results, which accord with the conclusions 
of Sizoo" as to remanence and induction but 
not as to coercive force, in general conform to 
the modified Weber theory which has been 
described by Stoner.” This theory, in brief, 
assumes that ferromagnetic substances contain a 
large number of domains. With fields greater 
than the coercive force but less than the satura- 
tion value each domain is assumed to have 
certain fixed directions which magnetization may 
follow in an attempt to approach parallelism 
with the impressed magnetic field. As the field 
is increased above saturation the directions are 
no longer fixed but may rotate so that parallelism 
to the field is then approached continuously. 


J. Sizoo, Properties of Ferromagnetic Crystals. 
Physica 10, 1-18 (1930). — 
2 E. C. Stoner, Magnetism, pp. 61-66 (London, 1930). 


uG, 


The curves for the two conditions of orienta- 
tion resemble those showing the ‘‘shearing effect” 
caused by a break in a solid substance," but in 
our tests the samples in both states of orientation 
contained the same number of joints with the 
same distance between particles (very nearly), 
since the number of particles and the packing 
density were the same. The effect of size of 
particles on magnetic properties in the size range 
used is small." 

If during orientation in a magnetic field the 
particles turn so that the directions of the re- 
sultants of their component directions of easy 
magnetization are parallel to the field, it would 
be natural that after demagnetization the first 
small increment of field strength in excess of the 
coercive force of the mineral should cause the 
greatest difference in induction per unit of field 
strength between the controlled and randomly 
oriented samples. The difference approaches zero 
as the field increases, as would be expected if 
the directions of easy magnetization become 
aligned in the field and the two samples become 
more nearly magnetically alike. The magnetiza- 
tion does not rise to a maximum at once or even 
when the coercive force is reached because the 
individual particles contain a large number of 


“domains” of magnetization, but if the mag- 


J. A. Ewing, Magnetic Induction in Iron and Other 
Metals, p. 279 (London, 1891). 
14 See reference 7. 
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netically oriented particles have a greater pro- 
portion of their ‘‘domains’’ in a position for easy 
magnetization than the random particles, the 
former would continue to show greater mag- 
netization than the latter for any given field 
until saturation was attained. 

Whether the slight irregularities of perme- 
ability exhibited by our data are due to the 
Barkhausen effect, to the movement of particles 
or to errors in the measurement of the current 
applied to the field magnet, has not been de- 
termined. 

On reduction of field strength the hysteresis 
curve for the aligned lies above that for the 
random particles, and the remanence in the 
former case is considerably higher than in 
the latter. One important factor contributing to 
this difference is apparently that the magnetized 
particles aligned with the north end of one 
adjacent to the south end of the next form 
elongated strings which reduce the demag- 
netizing effect as compared to the end effect of 
the short unordered particles. 

The fact that coercive force is independent of 
particle orientation indicates that the domain 
direction of easy magnetization is also a direction 
of easy demagnetization, for although the rema- 
nent magnetism may differ, the same opposed 
field strength effects demagnetization for differ- 
ent orientations. 

From the theory of Dean” and the results of 
Gottschalk'® the high coercive force of reduced 
hematite indicates a large interfacial area and 
the possibility of small ultimate particle size so 


15 See reference 3. 

‘°V. H. Gottschalk, The Coercive Force of Magnetite 
Powders, Report of Investigations, 3268, U. S. Bureau of 
Mines (in course of publication. 


that the smaller differences in magnetic prop- 
erties for different orientations of the reduced 
hematite may be due to a more random dis- 
tribution of its magnetic elements than in the 
crystalline magnetite. Such a condition is sub- 
stantiated by the work of Doan" who found by 
x-ray diffraction methods that for minus 200- 
mesh particles oriented in a magnetic field the 
magnetite in reduced hematite was randomly 
oriented while Mineville magnetite was not. 


CONCLUDING REMARKS 


Other work prevents an extension of the in- 
vestigation to determine whether the study of 
pure crystals by the method outlined may have 
more fundamental significance in the field of 
ferromagnetism. 

This method of studying the effect of particle 
orientation on magnetic properties offers certain 
advantages over those previously used. Some of 
the outstanding characteristics of the method 
are: 

(1) It permits the magnetic measurement of 
small samples of small particles, thus enabling 
production and use of very pure and uniform 
material. 

(2) It allows a single sample to be oriented in 
different ways without recourse to mechanical 
work or heat treatment, thus avoiding un- 
certain changes in chemical and structural com- 
position. 

(3) It permits the use of any size of particles 
desired without heat treatment and allows the 
direct microscopic measurement of the particles. 

The author wishes to thank Dr. V. H. Gott- 
schalk for his generous assistance and advice. 


1 —D. J. Doan, Private communication, 1934, 
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Comparison of the Temperatures in a Solid and Its Scaled Model 


R. V. CHURCHILL, Department of Mathematics, University of Michigan 


(Received December 1, 1934) 


If a solid is constructed of the same material and to 
scale with a body of any shape and subjected to initial 
and surface temperature conditions of the same character, 
properties can be derived to compare or contrast the 
temperatures at corresponding points in these solids. 
Seven properties are given in this paper, covering the cases 
of variable and steady temperatures of solids with and 
without radiation or convection at any part of the bound- 
ary, and with any fixed temperature distribution on the 
rest of the surface. The law used for the transfer of heat by 


1. INTRODUCTION 


F one body is a scaled reproduction of another 

the two are said to be geometrically similar. 
Two such solids can be placed in perspective 
(in an infinite number of ways) so that to each 
point (x,y,z) of the first there corresponds a 
point (ax, ay, az) of the second, when the origin 
is taken as the point of similitude. The number 
a is positive and constant for all points, and the 
second solid is a times as large as the first. 

When two such similar solids of the same 
material are subjected to the same initial and 
surface temperature conditions at corresponding 
points, their 
points generally have different values at any 


temperatures at corresponding 
instant. Comparisons of these temperatures are 
made in this paper. Considerable information is 
obtained about these temperatures without re- 
stricting the shape of the solid, even when the 
thermal coefficients are variable and the surface 
heat transfer law is quite general. A comparison 
of these temperatures in the case of steady 
temperatures in solids with constant thermal 
coefficients has already been carried out.' 

We consider here any two geometrically similar 
solids of arbitrary shape, made of the same 
homogeneous’ and isotropic material. The ther- 
mal conductivity can be any positive function 


1R. V. Churchill, General Properties of Steady Tempera- 
tures in Solids Partially Exposed to Gas, Physics 4, 50 
(1933). 

2See remarks on non-homogeneous solids in Section 5 
below. 


radiation and convection is a very general one. The thermal 
coefficients are functions of the temperature. These prop- 
erties are derived by substitutions of new variables into 
the heat equation and the boundary conditions. Some of 
the properties give important information about the effect 
of the size of a solid upon its rate of cooling or heating. 
All of the properties are useful in the determination of 
temperatures in any solid by measuring the temperatures 
in a large or small scale model, 


of the temperature, k = x(7); likewise the specific 
heat and density can be functions of 7. More- 
over, if there is radiation (or absorption) or con- 
vection of heat, or both, at any portion of the 
boundary surface the flux of heat through this 
portion at any point can be any function of the 
position and temperature of this point, 


Flux = ®(x, y, 2, T). (1) 


This portion of the boundary will be called the 
“exposed surface.’’ We define the flux as positive 
if heat flows into the solid through this surface. 

This surface heat transfer function ® is sub- 
jected only to the natural restriction that it 
cannot increase with 7 at any point. This is 
natural since it means that the absolute value 
of the flux cannot increase when the surface 
temperature is brought nearer to the temperature 
of the surroundings. In particular, this law of 
transfer includes Planck’s radiation formula and 
the Stefan-Boltzman law of radiation, as well as 
the law of Dulong and Petit for radiation and 
convection, provided the state of the surrounding 
gas and radiating or absorbing bodies is kept 
fixed. This law also includes the condition of 
insulation (x, y,s, 7)=0 as well as that of a 
constant supply of heat (x, y, sz, 7) =constant, 
over any part of the boundary. 

Any fixed temperature distribution may be 
assigned to the rest of the boundary surface, and 
the initial temperature distribution in the solid 
is arbitrary. 
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2. COMPARISON OF VARIABLE TEMPERATURES IN 
THE CASE OF RADIATION AND CONVECTION 


Consider a solid of any shape with an arbitrary 
initial temperature distribution. Let f(x, vy, s) =0 
be the equation of that portion of its boundary 
surface which has a fixed temperature at each 
point, and let g(x, y,s)=0 represent the re- 
mainder of its boundary at which heat transfer 
by convection, radiation or absorption, or a 
combination of these, takes place in accordance 
with the general law (1).* The temperature 7 at 
any point at any time is then determined by the 
conditions, 


c OT dt=(0 dx)(k OT Ax) 
+(0/dy)(k dT dy) (2) 
+(0/0z2)(k dT dz), 


T=T (x, y, z) when t=0, (3) 
T=T,(x, y, 2) on f(x, y, z) =0, (4) 
kdT/dn=9(x, y,2, 7) on g(x, y, 3) =0. (5) 


The heat capacity c per unit volume is any 
function of the temperature, and » is the out- 
ward drawn normal to the surface g(x, y, z) =0. 
Eq. (2) can be reduced‘ to the form of the 
heat equation with constant conductivity by 
the substitution of the new dependent variable 


u= | 
bela |) 


The conductivity «(7) is positive, so u is a 
monotonic increasing function of 7, and Eq. (6) 
defines T as a unique function of u, 


7 


K(A)dX. (6) 


T=r(u). (7) 


Since du/dT =x(T) =k, Eqs. (2)—(5) become, 
in terms of u, 


Ou 

r =(k/c)(0?u/dx?+0?u /dy?+0?u/d2), (8) 
c 

u=Uo(x, y, 2) when (=O, (9) 


u=u,;(x, y,2z) on f(x, y, 2) =9, (10) 


‘Either of the functions f and g may be represented by 
more than one analytic expression. For instance f(x, y, 2) 
=( may represent one base and the lateral surface of a 
cvlinder. 

* This has been pointed out by M.S. Van Dusen, Note 
on the theory of heat conduction, Bur. Standards J. Research 
4, 753 (1930). 
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du dn=(x, y, 5, r()) 


= (x, Vv, 5,4) on g(x, y,s)=0, (11) 
where 
in ¥. 2) 
Uo(xX, VY, 5) = K(A)dX, 
“0 
oTs(z, y, 2) 
u,(X,.¥, 3) = K(A)dX. 


Let the solution of this boundary value prob- 
lem be 


u= U(x, y, g, t). (12) 


We now substitute new independent variables 
throughout this problem, 


ol ce acta Faas t = ae 
x =ax, VY=ay, 2 =azZ, 


l’=a’t, (13) 


where a@ is a positive constant. If \, uw, v are the 
direction cosines of the normal », then 


du/dn=d0u Ox+pdu dy+vdu dz 


=a(dA\du dx’ +pdu ay’ +vdu/ds’). (14) 


We shall call the expression in parentheses 
du/dn'. 

When these new variables are substituted the 
factor a cancels in Eq. (8), and conditions (8) 
(11) become 


du’ (at =(k/c)(02u' Ax” +02u' dy” 


+0?u'/ds’) (15) 


uw’ =uy(x'/a, va,’ a) when U=0, (16) 


u’=u,(x’ /a, vy’ /a, 2’ /a) 


on f(x’ /a, v/a, 2’/a)=0, (17) 


du’ /dn' =(1/a)o(x' /a, y’ a, 2’ /a, u’) 


on g(x’ /a, y’/a,2’/a)=0. (18) 


Here w’ is the function of x’, y’, 2’, t’ into which u 
transforms under the substitutions (13), so 
according to Eq. (12) the solution of this prob- 
lem is 

u’ = U(x'/a, y’/a, 2'/a, t'/a?). (19) 


Now let x’, y’, 2’ be interpreted as the coordi- 
nates of a point, and ¢’ as time, in terms of the 
same units of length and time as those used 
originally in Eqs. (8)—(12). Eqs. (15)-(19) then 
have a new physical interpretation. 

A surface F(x/a, z/a)=0 is in per- 
spective with the surface F(x, y, z) =0, with the 


y/aQ, 
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constant a for its ratio of perspectivity, and with 
(x, y,2) and (x/a, y/a,2/a) as corresponding 
points. Hence Eqs. (17) and (18) are the sur- 
face boundary conditions for the function u in 
a solid in perspective with the original solid. 
A comparison of the conditions (17) and 
(10) show that the value of w’ at (x’, y’, 2’) on 
the surface f(x’/a, y'/a,2’/a)=0 is equal to 
uo(x’ /a, y'/a, 2’/a), which is the value of u at 
the corresponding point on the boundary of the 
original solid. It is easily seen that the direction 
cosines are the same at corresponding points of 
two surfaces in perspective. So, according to 
Eq. (14), du/dn’ is the normal derivative of the 
function u’ at the boundary of the new solid. 
A comparison of Eqs. (18) and (11) shows that 
the normal derivatives of u’ and u at corre- 
sponding points on the surfaces of the two solids 
vary with uw according to the same law except 
for the factor 1 /@ in Eq. (18). 

According to Eqs. (16) and (9) the initial 
values of u’ and u are the same at corresponding 
points of the two solids. Eqs. (15) and (8) show 
that u’ and u satisfy the same differential equation 
throughout the solids. Finally, a comparison of 
Eqs. (19) and (12) shows that the value of wu’ at 
any point of the new solid at time ¢ is equal to 
the value of u at the corresponding point of the 
original at time ¢ a’, 


‘la, t/a’). 


u’ = U'(x’, y’, 2’, ') = U(x'/a, y'/a, 2 

Since k and ¢ are the same functions of T in 
both solids, it follows from Eq. (7) that the 7”’s 
are equal when the w’s are equal; also if the u’s 
are unequal the 7’s are unequal in the same 
sense. The first property is thus established. 

Property I. If a solid 1s constructed to the scale a 
-with another and subjected to the same initial tem- 
- perature and surface temperature conditions at 
corresponding points, except that its surface heat 
transfer function is 1/a times the function for the 
original, the temperature at each of its points at 
time t is the same as the temperature at the corre- 
sponding point of the original at time t/a’. 

A less involved property arises from this if we 
assume that the thermal coefficients vary so 
that the ratio k/c is constant. This condition is 
closely approximated in various applications. 
Under this assumption a change in the character 
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of the surface heat transfer which causes the 
function ¢ to be replaced by ¢’, where 


¢' (x, ¥, 2, u) SG(x, y, 2, u), (20) 


will produce a new uw less than or equal to the 
original. 

To prove this we first write the conditions 
which determine u’, the new u: 


du’ /dt=(k/c)(@u' /dx*+6u' /dy" 


+6?u'/ds*), (21) 


(22) 
(23) 
(24) 


Conditions (8)—(11) determine u in the same 
solid. By subtracting each equation in the set 
(21)-(24) from the corresponding one in the 
set (8)—-(11) and writing w=u—v’, we find that 
w must satisfy the conditions 


u’=uo(x, y, 2) when ¢=0, 


u’=u;(x, y, 2) on f(x, y, z) =0, 


du’ /dn= 9'(x, y, 2, u’) on g(x, y, z) =0. 


dw /dt=(k/c)(Pw/dx?+0w/dy?+ Hw) ds"), 
w=0 when t=0, 
w=0 on f(x, y, z)=0, 


dw/dn= (x, y, 2, u) — (x, ¥, 2, u’) 


on g(x, y,z)=0. (25) 


These conditions show that w can be in- 
terpreted as the variable temperature in the 
given solid with initial temperature zero and 
temperature of the surface f(x, y, z)=0 held at 
zero, while heat transfer takes place at the 
surface g(x, y, 3)=0. However, if w=0 at any 
point of this surface then u=w’ there and 


o(x, vy, 3, w= O(x, y, 2, wu’), 


since @ is a non-increasing function of 7 and 
therefore of u. Hence 


o(x, y, 2, u)—o' (x, y, 2, u’) 
= o(x, y, z, u’)—o'(x, y, 2, uw’) 20, 


according to the condition (20). It follows from 
Eq. (25) that, at each point of g(x, v, s) =0, 


(26) 


dw/dn=0 when w=0. 


Let us now resort to physical argument. Since 
n is the outward drawn normal to g(x, y, 2) =0 
the condition (26) states that heat cannot flow 


out of the solid at any point on this surface 
whose temperature w is zero or lower. Since the 





TEMPERATURES 


initial temperature is zero throughout and the 
remaining surface f(x, y,z)=0 is held at zero 
temperature, it is evidently impossible for w to 
fall below zero at any point on the boundary. 
It follows from the second law of thermo- 
dynamics that w=0 throughout for all ¢. Hence 
when the condition (20) is satisfied, u’(x, y, 2, t) 
_ u(x, y, 2, t) throughout. 

We now apply this result to Property I. If 
d(x, ¥, 3, 4)=O on g(x, y, z2)=0, then (1/a) 
K o(x, ¥, 2, 4) SO(x, vy, 2,u) when a>1. Hence 
when the transfer function ¢/a obtains, the 
value of u is not greater than that resulting from 
the transfer ¢. A similar argument applies if 
o(x, ¥, 3, 4) =0, so we have proved 

Property II. If a solid is constructed to the 
scale a>1 with another and subjected to the same 
initial and surface temperature conditions at corre- 
sponding points, the temperature at each of its 
points at time t is greater than or equal to the tem- 
perature at the corresponding point of the original 
at time t/a? as long as the flux of heat through the 
exposed surface has not been outward at any point. 
If this flux has not been inward then the tem- 
perature is less than or equal to the corresponding 
temperature at t/a®. Here it is assumed that k/c 1s 
constant for these solids. 

A less accurate statement is the following: if 
the size of a solid is increased by a factor a, it 
heats more than 1/a® times as fast and it cools 
more than 1/a? times as fast. 


3. COMPARISON OF VARIABLE TEMPERATURES IN 
CasE THERE Is No RADIATION 
OR CONVECTION 


We now consider the important special case 
in which 
O(x,y,2,7)=0 on g(x, y,s)=0. (27) 
As before, k and ¢ are functions of 7, the initial 
temperature throughout the solid and the fixed 
temperature on the surface f(x, y,z)=0 are 
arbitrary, but the remainder of the boundary is 
insulated (condition (27)). Our next property 
follows directly from Property I. 
Property III. When part of the boundary surface 
is held at a fixed temperature distribution while the 
rest is insulated, an increase in the size of a solid 
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by a factor a causes an increase in its heating or 
cooling time by a factor a’. 

More precisely, if a solid is constructed to the 
scale a with a given solid and subjected to the 
same initial and surface temperature conditions 
at corresponding points, the temperature at each 
point at time ¢ equals the temperature at the 
corresponding point of the original at time t/a’. 
Thus it takes this solid a? times as long to attain 
a temperature at each point equal to that at the 
corresponding point of the given solid. 

In particular, if the temperature at any point 
in the given solid is always increasing, the 
temperature at the corresponding point in the 
solid @ times as large at the same instant is lower 
if a>1. 

Property III also applies to solids whose entire 
boundary is subjected to a fixed temperature 
distribution (condition (4)). 


4. COMPARISON OF STEADY TEMPERATURES 


The steady temperature is the limit of our 
variable temperature T, as to. If this exists 
at each point of the solid then the steady value 
ui of the function wu is obtained from Eq. (6), 


? 


K(A)dX. 


“0 


u=lim 


too 


By interpreting x’, y’, 2’ again as the coordi- 
nates of a point in a solid a times as large as the 
original, it follows from Eq. (19) that 


ai’ = lim U(x'/a, y’/a, 2’/a, t/a?) 


t'—>c0 


= lim U(x'/a, y’/a, 2'/a, t). 
t 


—coO 


According to Eq. (12) this second limit is the 
steady value of u at the corresponding point of 
the original solid, so #’=a at corresponding 
points. If ¢ is continuous then it follows that 
lim $(x’/a, y’/a, 2’/a, u’) 

t'—»00 


=lim $(x'/a, y’/a, 2'/a, u); 


t—>co 


so the steady values of the transfer functions are 
equal at corresponding points, except for the 
factor 1/a@ in condition (18). Hence for solids 
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thermal coefficients functions of 
temperature, we have 

Property IV. When the size of a solid is in- 
creased by the factor a with surface temperature 
conditions kept the same at corresponding points 
except that the transfer function is multiplied by 
the steady temperatures at corresponding 
points are equal. 

The steady function @ satisfies Laplace's 
equation, since da, dt=0 (Eq. (8)). But @ is a 
monotonic increasing function of the steady 
temperature, so the properties of steady temper- 


whose are 


i ay 


atures already derived! apply to @ and hence to 
steady temperatures in with variable 
thermal coefficients. We state the two of greatest 
importance below. 

Property V. An increase in the size of a solid 
by a factor a produces the same change in the steady 
temperature (at corresponding points) as that 
caused by multiplying its transfer function @ by a. 

That is, if a body is constructed of the same 


solids 


material and to the scale a@ with a given one, 
with the same fixed temperature at corresponding 
surface points and the same transfer law at all 
other corresponding points on the surface, its 
steady temperature at each point equals that 
which would be obtained at the corresponding 
point of the original by multiplying the function 
¢ of the original by a. 

Property VI. An increase in the size of a solid 
causes an increase in its steady temperature at 
corresponding points if the solid is absorbing heat 
at its exposed surface (@>G@), and a decrease if it 
is losing heat through this surface. The fixed 
fem perature surfaces are excepted, of course. 

In case there is no radiation or convection at 
the surface, (x, y,2,%7)=0 on 
Property IV then becomes 

Property VII. When two solids are constructed 
to scale with corresponding parts of their boundaries 
held at the same fixed temperatures and the rest 
insulated, the steady temperatures at corres ponding 
points throughout are the same. 


g(x, v, s) =0. 


5. Discussion 


These properties appear to be useful in a 


variety of problems of design involving large or 
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small scale models when it is desirable to 
determine temperatures from these models. 

Property I, the most general of the above 
properties, enables one to determine the temper- 
ature at any point in a solid at any time from 
the temperatures in a scaled model, provided 
the transfer of heat at the exposed surface can 
be properly controlled. If this control is not 
possible so that the same law of heat transfer at 
the exposed surface must be applied to the model 
and the original solid, then Property II can be 
used. This furnishes a lower limit for the temper- 
ature at each point of a solid in terms of a 
corresponding temperature in a small scale 
model, provided the solids are absorbing heat 
through their exposed surfaces. It furnishes an 
upper limit if the solids lose heat through these 
surfaces. 

The application of the rest of the properties, 
dealing with insulated surfaces and_ steady 
temperatures, to determine the temperatures in 
a body from those in a model is more direct and 
quite obvious. 

The condition (19) or its equivalent, 


U'(x, y, 3, t) = U(x/a, y/a, 2/a, t/a*), 


shows that the u’=constant surface at time ¢ a? 
is in perspective with the “=constant surface 
at time ¢, with @ as the ratio of perspectivity. 
Hence in a scaled model satisfying the conditions 
of Property I or Property III, every isotherm in 
the model corresponds to a geometrically similar 
one in the original solid; but the temperatures 
on these isotherms are equal at different times. 

Under the conditions of Property IV or Prop- 
erty VII (steady temperatures) each isotherm 
in the model is geometrically similar to the 
corresponding isotherm in the original solid. 
The ratio of magnification is the same for the 
isotherms as for the solids. 

For composite bodies consisting of two or more 
homogeneous solids in perfect thermal contact, 
Property I still holds true under the same general 
initial and boundary conditions. Property III is 
also true for composite solids since it is a special 
case of Property I. The proof will be given ina 
later paper in which a more complete investiga- 
tion of the problem for non-homogeneous solids 
will be made. 
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rolled and in cold rolled sheets of silicon steel 

ire directions Of easiest magnetization both before 

ifter the final high temperature anneal. They are 
exhibited clearly by a disk suspended in a magnetic field. 
This property can be explained by the fact that a great 
part. of the crystals composing the sample lie in a preferred 
orientation. The samples thus exhibit to a higher or smaller 


degree the anisotropic magnetic properties of single 
crystals. This non-random orientation of the crystals was 
shown for the annealed samples by a determination of the 
orientation of a great number of crvstals bv a light re- 
flection method. The preferred orientations in cold rolled 


and hot rolled annealed samples differ substantially. 





GENERAL SURVEY 


ARLY in the development of silicon steel as 
used in transformers and rotating electrical 
machinery, it was recognized that the sheets 
showed a nonuniformity in their magnetic quali- 
ties. With one and the same sheet, it was found 
that permeability and hysteresis loss had differ- 
ent values depending on the direction in the 
plane of a sheet in which they were measured. 
This was true both immediately after the sheets 
had been rolled to their final thickness and also 
after the following final anneal which improved 
the magnetic properties in general. The non- 
uniformity resulted in the requirement in the 
Epstein test for magnetic losses that half of the 
total material tested must be measured with the 
flux in the rolling direction, and half of it with 
the flux perpendicular to this direction. Thereby 
a certain average of the qualities of the material 
was obtained. 

Recently. Dahl and Pfafienberger' have em- 
ployed a very simple method by which the 
inhomogeneous properties of magnetic samples 
can be demonstrated in a convincing way. A 
disk is punched from the sheet which is to be 
investigated and suspended in a magnetic field 
of sufficient strength. If this disk were homo- 
geneous (isotropic) magnetically, it would re- 
main in any position in the field. If, however, 
the material has a higher permeability in one 
direction than in another, the disk will assume 
a position in which this direction of high per- 
meability is parallel to the field lines. In any 
other position a torque will be developed which 


'E. Dahl and F. Ptaffenberger, Zeits. f. Physik 71, 93 
1931). 


tends to pull 
positic m. 


the disk into its equilibrium 


We are here mainly concerned with the cause 
for this inhomogeneity, which, due to the ad- 
vance in our knowledge of magnetism during the 
last 10 years, seems now to be definitely estab- 
lished. Several explanations can be and at various 
times have been advanced. The fact that the 
grains are elongated in the rolling direction 
might explain a smaller loss in that direction 
than perpendicular to it because there the flux 
has to pass fewer obstructing grain boundaries. 
It would, however, hardly explain the fact that 
in certain cases the losses are lowest at 45° to 
the rolling direction, and it certainly would not 
hold for a sheet after anneal where the grain 
dimensions are, on the average, the same in all 
directions. The same arguments rule out internal 
stresses as the cause for the nonuniformity, 
although they will probably have some effect 
before annealing. Messkin and Kussmann? as- 
cribe the directional properties to non-metallic 
inclusions elongated in the rolling direction by 
the process of rolling which would not be affected 
by the annealing. But this effect, just as the 
first one mentioned, could not explain the 
presence of a lowest loss direction at 45° to the 
rolling direction. 

The correct interpretation was made possible 
through our advanced familiarity with the mag- 
netic properties of individual crystals of ferro- 
magnetic materials. Little was known about 
them as long as one dealt oniy with aggregates 
of many small crystals such as are met with in 
commercial steel. Later, by special methods, it 


2W.S. Messkin and A. Kussmann, Die Ferromagnetischen 
Legierungen, 1932. 


105 








106 K. Jj. 


was possible to produce large crystals on which 
magnetic measurements could be _ performed. 
Ruder,’ in 1917, made strips of Si steel of up to 
10 inches in length, which in their entire length 
consisted of one single grain. Such single crystals, 
as they are commonly called, have different 
magnetic properties in different directions just 
as rolled sheets have, and a possible connection 
suggested itself. If it could be proved that the 
small crystals in rolled material were not oriented 
at random, but showed certain preferred orienta- 
tions, then the observed anisotropy could easily 
be explained as being caused by the magnetic 
anisotropy of the individual crystal. This proof 
was furnished for unannealed material, where we 
deal with a fine grained fiber structure, by an 
x-ray analysis using the Debye-Scherrer-Hull 
method. In the case of annealed material, how- 
ever, where the grain size is too large, x-ray 
analysis can, in general, give us no information 
as to whether or not a preferred crystal orienta- 
tion is present. Goss‘ has recently published 
x-ray Laue diagrams of annealed cold rolled 
specimens, and has concluded without justifica- 
tion, in the light of our results given below, that 
no preferred orientation was present. An optical 
method in which the crystal orientation is 
determined by light reflections from the etched 
surface, and which will be described later, showed 
quite definitely that a preferred orientation was 
also present after anneal in hot and cold rolled 
sheets. 


MAGNETIC PROPERTIES OF SINGLE CRYSTALS 
OF IRON 


At this point we stop to consider the magnetic 
properties of single crystals, firstly to show some 
properties of the crystalline state, and secondly 
to analyze what information could be obtained 
from measurements with the simple torque 
apparatus mentioned in the beginning. 

Instead of Si steel we choose pure iron whose 
properties are better known, but which, for our 
purpose, is representative of the former. In an 
iron crystal the iron atoms have a regular 
distribution—they form a so-called body cen- 


*W. E. Ruder, Trans. Am. Soc. Steel Treating 8, 23 
(1925). 

*N. P. Goss, Trans. Am. Soc. 
Meeting, 1934. 
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Fic. 1. The magnetization curves of iron along three 
crystal axes which are indicated on the elementary cube 
shown as insert. 


tered cubic lattice. The fundamental element of 
this lattice is a cube whose corners and whose 
center are each occupied by an iron atom. (Fig. 
1.) It was found by measurement that up to 
about J = 1000 (B~ 12,000 gauss) the magnetiza- 
tion curve practically coincided for all directions 
in the crystal, but that above that value the 
magnetization was different depending on the 
direction in which H was applied. It took only 
small fields to saturate the crystal along a cube 
edge ({ 100] axis), and increasingly higher // for 
saturation along the face diagonal ([110] axis) 
and the body diagonal ({111 ] axis). Accordingly 
the six [100] axes are called the axes of easiest 
magnetization. According to our present view 
the process of magnetization along a [110] axis, 
for example, takes place in the following manner. 
The magnetization vectors in the various ele- 
mentary districts of which the crystal is built 
up gradually, upon increase of H, grow in the 
two [100 ] directions closest to the field direction 
until half of all the districts are magnetized in 
each of these two directions, and the component 
of I in field direction is J,/\2. This process 
requires practically no energy, but with further 
increase of // the I vectors swing gradually into 
the field direction, and the energy required is 
much greater. This explains the bend in the 
magnetization curve for the [110] axis at 
I=I,/\2. The magnetization along the [111] 
axis takes place in a similar way, the bend 
appearing at J=J,/ 3. The so-called technical 
saturation J, is reached when I lies parallel to //. 





$ 
. 
y 
» 
r 
) 
y 
t 


MAGNETIC ANISOTROPY IN SILICON STEEL 107 


It is quite evident from these facts that a 
single crystal disk, having different permeabilities 
in the different crystal directions, would orient 
itself in a magnetic field which is sufficiently 
strong to produce an induction above the first 
bend. The same follows logically for a polycrystal 
whose crystals are not oriented at random, but 
show certain preferred orientations. The mathe- 
matical expression for this torque can be derived 
easily under certain simple conditions. 

Akulov® first set up the fundamental energy 
equation for ferromagnetic crystals, free of 
strains, expressing their free energy F in terms 
of the direction which magnetization I and field 
H have in the lattice. 


F= Fy(1) +2k(en?a2? + a22a3? + a132011") 
— HI (0181+ a282+ 383). (1) 


F,(I) is the energy necessary to magnetize the 
crystal along a 100 axis up to the intensity J, 
kis aconstant, a, a2, a3 are the direction cosines 
of the magnetization I with respect to the three 
main crystal axes (J,/I,, -++) and #1, Be, 83 the 
direction cosines of H with these axes (//,/H, 
-++). Because of the second term on the right 
side of the equation the direction of I is in general 
not the same as that of J/. In a disk, cut from 
this crystal and suspended in a field, this devia- 
tion will result in a torque, and if we take the 
particular case of a disk cut parallel to a (100) 
plane, we can derive the torque T from Eq. (1) 
by differentiation : 


dF/dg=T=k sin 4¢( =HI,). (2) 


y is the angle between magnetization and a [100 ] 
crystal axis, J, the component of I normal to HZ. 
The constant k may be derived by measuring 
H and I, or T.* We can easily see that in a 
polycrystal, a certain percentage of whose crys- 
tals had a preferred orientation with a (100) 
plane in the surface and the rest of whose 
crystals were distributed at random, we would 
obtain a torque curve of the type of Eq. (2), 
only with a different constant k’. k’/k would 
then represent the percentage of crystals in 
the preferred orientation. N. Akulov and N. 


5 N. Akulov, Zeits. f. Physik 69, 78 (1931). ; . 
* The torque due to rotational hysteresis is negligible 
at the high field strengths used. 


Briichatov’ have used this method for the 
analysis of rolled sheets and derive the percentage 
of crystals in certain preferred orientations 
without, however, offering a check by some other 
method. In general this way of analyzing is only 
possible if the direction or directions of preferred 
orientation are known from the start, and if the 
spread about this orientation is small. Thus it 
becomes necessary, before a quantitative inter- 
pretation of the torque curves can be under- 
taken, to establish the direction and degree of 
preferred orientation by an x-ray method, or, if 
that is not possible, by a light reflection method. 
Although a quantitative interpretation of the 
torque measurements thus necessitates the taking 
of additional data, the torque method in itself 
is an excellent and extremely simple means of 
obtaining qualitative information about the 
magnetic merits of a material. The greater the 
torque maximum, the larger will be the number 
of crystals that have cube edges ([100] axes) in 
a uniform direction. The more crystals with cube 
edges there are in this direction, the higher will 
be the permeability, again in this direction, at 
high fields. (See Fig. 1.) The angle between this 
preferred direction and the rolling direction can 
be read from the torque curve as the equilibrium 
position next to the torque maximum under 
consideration. 

It seems advisable at this place to discuss 
briefly the advantages which a material with 
preferred orientation has over materials with 
random orientation, especially since a fine grain 
material with uniform properties in all directions 
was in the past considered as a condition which, 
although it could not be achieved, should at 
least be approached. The most important prop- 
erties of silicon steel as used, for instance, for 
transformers are its permeability and its hys- 
teresis loss. For the former one can make 
definite predictions on the basis of Fig. 1. 
Below a magnetization of J~1000 (B~12,000 
gauss) the permeability is practically the same 
in all crystal directions and only little can be 
gained by having a preferred orientation. Above 
this value, for B>12,000 gauss, the increase in 
permeability through preferred orientation—if 


7N. Akulov and N. Briichatov, Ann. d. Physik 15, 741 
(1932°. 
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the field is always in the direction of this pre- 
ferred orientation—is very marked and a given 
B-value can be reached with much smaller // 
than without a preferred orientation. Little can 
be said about the effect of preferred orientation 
on hysteresis. In a series of measurements on 
single crystals a dependence of hysteresis on 
crystal orientation was found® with a minimum 
along the [100] axis, which would indicate that 
a preferred orientation might result in a reduc- 
tion of hysteresis. The effect of preferred orienta- 
tion on eddy current losses can reasonably be 
expected to be very slight. 


TORQUE MEASUREMENTS 


Fig. 2 shows the apparatus developed for 
determining the torque at any angle between 
field and rolling directien. A scale is attached to 
the disk to be measured and when the movable 
pointer P is turned, it will indicate the torque 
in degrees at any position of the disk which in 
turn is observed at the stationary pointer A. 


e tera 


Electromagnet 


Fic. 2. Apparatus for measuring the torque exerted on 
a disk in a magnetic field (apparatus in an early stage of 
development). 


By calibrating the spring the torque can be 
expressed in gem or ergs instead of in degrees. 
The field value chosen had such a value that on 
further increase of the field the torque curve 
remained the same. Fig. 3 shows results obtained 
in two disks with different treatment. In a hot 
rolled sheet’ there are two preferred directions 
at approximately +45° to the rolling direction 
both before and after annealing. In a cold rolled 


*W. E. Ruder, Trans. Am. Soc. for Metals 22, 1120 
(1934). 

® The silicon steel used in all these experiments contained 
about 3.5 percent Si and 0.05 percent C. 
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Fic. 3. Torque curves at high fields measured with the 
apparatus shown in Fig. 2. 3a gives the curves for a hot 
rolled sheet, 3b for a cold rolled sheet, both before and 
after anneal. H =900 oersteds. 


sheet with 85 percent cold reduction,’ on the 
other hand, a distinct change takes place at the 
anneal. Before anneal the preferred direction is 
similar to the one found in hot rolled sheets, but 
after the anneal at 1100°C the recrystallized 
sheets have preferred directions which are paral- 
lel and perpendicular to the rolling direction. 

This result, which shows the typical difference 
between hot and cold rolling found in practically 
all the samples investigated, could be expected 
from x-ray measurements as far as the pre- 
anneal sheets were concerned. In these sheets 
x-rays reveal a preferred crystal orientation with 
a [100] axis at 45° with respect to the rolling 
direction. For the annealed samples a_ light 
reflection method was used to indicate the 
preferred orientation. 


LiGHTtT REFLECTION METHOD 


The sample surface is polished very carefully 
and etched for about 20 minutes in a 3 percent 
solution of ammonium persulfate. If the sample 
surface after‘etching is viewed in parallel light, 


10 See U. S. Patent 1915766, Smith et ai. 





MAGNETIC ANISOTROPY 


Fic. 4. Etched surface of a sample of cold rolled Si steel 
(20 magnification). 4a, vertical illumination; 4b, the 
same area with horizontal illumination with the light beam 
perpendicular to the rolling direction, 


a number of crystals appear bright while others 
are dull. This is due to the peculiar way in which 
the crystals are attacked by the etching agent. 
The latter opens its attack on a crystal in a 
great number of places digging pits into the 
surface! whose walls (and this is the important 
feature) are nearly plane and parallel to cube 
faces, i.e., (100) planes. If we illuminate an 
etched sample vertically under the microscope 
all the crystals which have (100) planes parallel 
to the surface accordingly will appear bright 
(=group I), the rest dull (Fig. 4a). Group I 
appears dark with horizontal illumination and 
stays dark as the sample is rotated through 360°. 
The rest (=group II), except a few that stay 
dark with both vertical and horizontal illumina- 
tion, go through a cycle of brightness followed by 
darkness during a 360° rotation, most of them 
twice, few of them only once (Fig. 4b). The 
group II crystals have very nearly (110) plane 
parallel to the sample surface, and, at the posi- 


"G. Tammann, Lehrbuch der Metallkunde, p. 129 (1932). 


SILICON STEEL 


—rolling direction 


Fic. 5. Three crystals on the same sample as shown in 
Fig. 4 under higher magnification (170). 


tion of maximum brightness, have a [100] axis 
very nearly perpendicular to the incident light 
if the two brightness maxima are close to 180° 
apart from each other. 

We still do not know the orientation of the 
crystals with a (100) plane in the surface. This 
can be found with a microscope using a high 
magnification with which the orientation of the 
individual etching pits may be observed. These 
pits, Fig. 5, form almost perfect squares on which 
the direction of a [100] axis can be directly 
measured. 

Any preferred crystal orientation can thus be 
discovered by determining the orientation of 
each individual crystal among a sufficiently large 
number on the specimen provided the crystals 
are fairly uniform in size. Such counts covering 
an area with about 250 crystals on each sample 
have been carried out by Mr. W. C. Johnson in 
this laboratory, and two representative deter- 
minations are plotted in Figs. 6 and 7 referring to 
the same samples which were used for Fig. 3. 
Incidentally, these counts have not only been 
made on the surface, but also in inner layers of 
the sheet after etching it to about half thickness. 
There, with slight variations, the same distribu- 
tion curve was found. All the crystals referred 
to as (100) crystals reflected light vertically and 
had a (100) plane parallel to the surface with a 
spread of about +10°. All crystals reflecting 
nearly horizontal light into the microscope are 





&: 2B; 


(110) crystals 
(55.3%) | 
> 
ee 


S 
é 


Sang 


ge 


4 





Percenta 


Random 
orientation: 19.6% 


Fic. 6. The distribution of [100] axes in (100) crystals 
(right) and in (110) crystals (left) on the hot rolled, 
annealed sample of Fig. 3a. The radius of each 10° sector 
gives the percentage of crystals having a [100] axis 
within these 10°. The inserts show the preferred orienta- 
tions. 
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Fic. 7. The distribution of [100] axes in (100) and (110) 
crystals on the cold rolled annealed sample of Fig. 3b. 


classified as (110) crystals if they showed two 
reflections at nearly 180° on rotation through 
360°. The deviation in this case, however, is 
greater, due to greater irregularities in the 
etching pit$, and the (110) plane in a (110) 
crystal may actually make an angle of up to 
+20° with the sample surface. This procedure 
of classifying together crystals of rather diverse 
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orientations is justifiable in our case since all of 
them have a [100] axis in or nearly in the 
surface, which is the essential property as regards 
magnetic behavior. All the crystals with only 
one and those with two reflection maxima which 
enclosed an angle <160° during a 360° rotation 
and the crystals remaining dark with both 
vertical and horizontal illumination are counted 
in the “random” class. In the hot rolled, annealed 
sample the (100) crystals have a decided pre- 
ferred orientation with a [100] axis at +42° to 
the rolling direction, whereas the (110) crystals 
show a slightly preferred orientation with the 
[100] axis in the rolling direction. In the cold 
rolled and annealed sample both the (100) and 
(110) crystals prefer a [100] axis in rolling 
direction. This sample has a particularly great 
uniformity in crystal orientation on account of 
an intermediate anneal.‘ 


DISCUSSION OF RESULTS 


These results are significant in two respects. 
From a metallurgical viewpoint they prove 
definitely that in iron with 3.5 percent silicon, 
after a high temperature annealing, the crystals 
are not oriented at random but show a preferred 
orientation whose direction and degree of per- 
fection depend on the preceding treatment. Hot 
rolling and cold rolling thus produce a strikingly 
different effect. Although their preferred orienta- 
tion is nearly the same immediately after the 
rolling process, the orientation in the cold rolled 
sheet—probably due to effects on slip planes and 
to strains—is entirely different from the hot 
rolled sample after recrystallization. Secondly, 
from a magnetic viewpoint, the findings agree 
qualitatively with the torque measurements and 
they also explain readily some of the magnetic 
properties. In hot rolled sheets, e.g., among the 
three directions at 0°, 45° and 90° with respect 
to the rolling direction, the highest permeability 
is found at 45°, with the 0° and 90° direction 
as second and third best. In cold rolled sheets 
the 0° direction is superior to the 90° and the 
45° one, and in both cases the preferred orienta- 
tion among the (100) crystals is the cause. 

Our findings enable us to test the accuracy of 
the method proposed by Akulov and Briichatov’ 
for deriving the percentage of crystals in the 
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preferred orientations from a torque curve. This 
was done by using the equations given in their 
paper which are derived from the fundamental 
Eq. (1), with the constant k=1.5 10° ergs cm~*. 
This value was determined by torque measure- 
ments on a single crystal disk of silicon steel. 
The results are given in Table I with the per- 


TABLE I. 








In surface _ (100) plane (110) plane 
In rolling direction [100] axis T110] axis [100] axis 


16(21) 
none 





Hot rolled sample 
Cold rolled sample 


none 
0(40) 


29(23) 
66(35) 








centages derived from our crystal count added 
in parentheses. 

The agreement for the hot rolled specimen is 
satisfactory, whereas no check is obtained in the 
case of the cold rolled sample. The latter result 
shows clearly that two different crystal distribu- 
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tions may result in the same torque curve. It is 
difficult to understand this possibility at first. 
If, however, one takes account of the possible 
spread of the (100) and (110) planes given above, 
the two preferred orientations given in Fig. 7 
are very nearly equivalent to one with a [100] 
axis in rolling direction and random distribution 
of the other [100] axes, and the result thus 
becomes more plausible. The two examples show 
that a quantitative interpretation of the torque 
curves is not possible in general. 

It was the main object of this paper to 
describe the two methods used for the investi- 
gation of magnetic directional properties. In a 
following article quantitative results on ma- 
terials after different treatments will be given. 

In conclusion the author wants to express his 
thanks to Mr. W. E. Ruder for many discussions 
on the subject and to Mr. W. C. Johnson for 
making the tests. 
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Growth of Crystals of Zinc Containing Cadmium by the Czochralski-Gomperz Method 


Haro_p K. SCHILLING, Physical Laboratory, University of Iowa 
(Received November 19, 1934) 


In previous experiments it was found impossible to 
duplicate results of Hoyem and Tyndall who found a 
“range”’ for the successful growth of single crystals of zinc. 
Most specimens were ‘‘optically mosaic.” It is now shown 
that by adding Cd to Zn of highest purity, single crystals 
may be grown and that a “range” is connected with their 
growth. The range is similar to, but not identical with, 
that of Hoyem and Tyndall. Moreover, the shape of the 
range is dependent on the Cd content and on the speed of 


INTRODUCTION 


NVESTIGATIONS of the conditions govern- 
ing the growth of zinc single crystals have 
been made in this laboratory by Hoyem 
and Tyndall,! Schilling? and Cinnamon.* The 
Czochralski-Gomperz method was used for the 


oi G. Hoyem and E. P. T. Tyndall, Phys. Rev. 33, 81 
1929). 


*H. K. Schilling, Physics 5, 1 (1934). 

*C. A. Cinnamon, R. S. I. 5, 187 (1934); W. J. Poppy 
has recently checked and extended Cinnamon’ s results, 
Phys. Rev. 46, 815 (1934). 


growth of the crystals. Optical mosaics as well as the 
changes of Hoyem and Tyndall are found outside of the 
range. In general optical mosaics are less frequent the 
higher the Cd content. Certain surface markings appear on 
the crystals even for only 0.20 percent Cd. These are 
believed to be identical with markings previously described 
by Rosbaud and Schmid and by Straumanis and probably 
indicate concentration of Cd along basal planes. 


first two investigations and a modified Kapitza 
method for the last. Hoyem and Tyndall found 
that single crystals of zinc, 2.7 mm in diameter, 
10 cm long, whose orientations were predeter- 
mined by ‘“‘inoculation’’ with a nucleus, could 
be grown at the speed of 1.2 cm/min. only under 
very restricted temperature conditions. It was 
necessary for them to keep the temperature of 
the molten zinc between upper and lower limits, 
only 4°C apart for orientations between 0° and 
55° but widening out rapidly for higher orienta- 
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tions. Outside this region’ of successful growth 
a crystal which started successfully’ changed to 
one of different orientation. They concluded that 
the determining factor in the successful growth 
of a crystal of a given orientation was a suitable 
temperature gradient in the liquid column below 
the growing crystal. 

Later Hoyem grew within the H—-T range a 
large number of successful crystals of ‘‘Kahl- 
baum” zinc.® He also found that outside the 
range failures (change to a crystal of different 
orientation) resulted. Immediately afterwards 
he also grew crystals of ‘‘spectroscopically pure”’ 
zine within the same range.’ It seemed desirable 
to investigate this matter further for other speeds 
of growth, sizes of crystals, etc. From the 
“successful ranges’’ which were expected to be 
found it was hoped to throw more light on the 
process of crystal formation. In attempting this, 
however, the writer found himself quite unable 
to reproduce the H—T results. Thinking that 
differences in the nature and amounts of im- 
purities in various brands of zinc might well 
influence the conditions of crystal growth, alto- 
gether nine brands and lots of zinc ranging from 
C.P. grade to grades certainly less than 0.01 
percent impure,* were used—all with the same 
results. Instead of single crystals optical mosaics” 
were obtained mainly. There were no ‘‘sudden”’ 
or “gradual’’ changes and therefore no range 
whatever. No changes in technique, nor use of 
other speeds of growth? altered the general 
character of the results. Further, two independ- 
ent observers, Professor E. P. T. Tyndall and 


‘ This will hereafter be called the H—T range. 

*Crystals were considered to have started growth 
successfully if they grew single and with the predetermined 
orientation for at least 1 cm. The chief enemy to a success- 
ful start was the failure to get a good junction between 
the nucleus and the new crvstal due to inclusion of oxide. 

* A. G. Hoyem, Phys. Rev. 38, 1358 (1931). That these 
crystals were certainly good has been shown by later 
evidence. See Poppy, reference 3. 

* These were considered single at the time but were 
under some suspicion on account of their anomalous 
electrical resistivities. See E. P. T. Tyndall and A. G. 
Hoyem, Phys. Rev. 38, 820 (1931). Lately reexamination 
has shown that some of these certainly were mosaic. See 
Poppy, reference 3. These S. P. crystals and the Kahlkaum 
crystals immediately preceding them were grown in a 
non-oxidizing atmosphere, a fact which must be taken 
into account in making comparison between their growth 
and the growth of crystals in air, as in the present work. 
_ ‘Ina previous publication, reference 2, the maximum 
impurity was stated erroneously to be 0.001 percent. The 
above statement gives the correct figure. 
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Dr. A. W. Hanson, kindly checked this conclu- 
sion using some of the writer’s material. 

Meanwhile, Poppy* and Hanson,’ using a 
different method, had obtained definitely non- 
mosaic crystals of lots of zinc with slightly less 
than 0.01 percent impurity, and Cinnamon* had 
even established for successful growth of single 
crystals of the same sort of zinc, conditions very 
similar to the H-T range. All these workers 
found, however, that optical mosaics were likely 
to occur for very pure zinc and that special 
precautions* were then necessary to prevent 
mosaic starts before single crystals could be 
grown. 

As the writer’s zincs were not all of high 
purity, it was not possible to conclude that 
optical mosaics were characteristic solely of very 
pure zinc, nor could it be supposed that the H-T 
range and the production of single crystals were 
simply characteristic of impure zinc, because it 
should then have appeared for some of the 
writer’s zincs. It seemed necessary, rather, to 
suspect it of being characteristic of some one 
particular impurity which might have been 
absent from all the writer’s zincs but present in 
the H—T zincs and later in the Kahlbaum zinc 
of Hoyem. Comparison by spectroscopic analy- 
sis!” showed definitely that all the zincs which 
had by the C-—G method yielded mosaics, had 
less (in many cases far less) Cd impurity than 
the lots for which the H—T range was established 
and the Kahlbaum zinc later used by Hovem. 

This led to the hypothesis that the mosaics 
might be characteristic of zinc containing no (or 
very little) Cd, but that by adding Cd to such 
material single crystals might be grown by the 
C-G method and that the H-T range might 
even be duplicated for some particular percentage 
of Cd. Preliminary experiments by Professor 
Tyndall showed that by adding two or three- 
tenths of one percent of Cd to very pure zinc"! 
the mosaics were indeed largely prevented, that 
the ‘“‘sudden”’ changes returned and that signs of 
a range of successful growth again appeared. 
This paper reports subsequent investigation of 
the matter for various additions of Cd to Zn. 


* A. W. Hanson, Phys. Rev. 44, 329 (1933). 

! Carried out by Professor E. P. T. Tyndall. 

" Horsehead-Special zinc produced by the New Jersey 
Zinc Company; also used by the writer in all of the present 
work. 
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APPARATUS AND PROCEDURE 


The apparatus was the same as that used 
during the writer’s earlier work.2 The growing 
procedure was essentially that of Hoyem and 
Tyndall. However, instead of using a stream of 
CO, to keep the surface of the molten zinc free 
of oxide, a speck of solid ZnCl,."* was dropped 
onto the surface just before inserting the nucleus. 
This melts and dissolves the surface oxide, 
assuring a good junction between the nucleus 
and new crystal. The deliberate forming of a 
neck at the start of a crystal was found un- 
necessary, and the writer began drawing out the 
crystal as soon as the desired temperature of 
growth was reached. The crystals were between 
2.5 and 3.3 mm in diameter and at least 10 cm 
long when grown at the speed of 1.2 cm/min. 
and 8 cm long at the speed of 0.55 cm/min. 


RESULTS AND DISCUSSION 


The results obtained with 0.26 percent Cd in 
zinc at the speed of growth of 1.2 cm/min. (the 
H-T speed) are shown in Fig. 1, abscissas being 
orientations and ordinates t—t, where ¢ is the 
temperature in degrees C of the zine directly 
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Fic. 1. Range obtained with 0.26 percent Cd+Zn at 
speed of 1.2 cm/min. @ single crystals; © failures by 
change from single to single of different orientation; + 
mosaics with no change of orientation; X mosaics with 
change of orientation. 





























" This was first suggested by W. Fahrenhorst (see E- 
Schmid and M. A. Valouch, Zeits. f. Physik 75, 532 (1932), 
who, however, used ZnCl, solution. 

See reference 1, Fig. 1. 
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TABLE I. Classification of specimens of Fig. 1. 





Syvm- 

bol 
in Near 
Fig. junc- 
Group 1 Type _ tion 
A @ a single 
b mosaic 
¢ mosaic 


Below 
junction 
single 
single 
improving 
mosaic 


Later change 


none 
none 
becomes single 


single 
mosaic 


discontinuous 
discontinuous 


single 
single 


single 
mosaic 


single 
single 


slowly contin. 
slowly contin. 


single 
mosaic 


single 
single 


rapidly contin. 
rapidly contin. 


mosaic 
mosaic 


single 
mosaic 


none 
none 


discontinuous 
discontinuous 


mosaic 
mosaic 


single 
mosaic 


mosaic 
mosaic 


single 
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below the liquid column under the growing 
crystal and f is the melting point of zinc. The 
black circles refer to single crystals. ‘‘Boundary”’ 
curves have been drawn in such a way that all 
these points, with one exception, are either 
definitely inside them or within 1° of them." 
The open circles, crosses and X's refer to failures 
of three types to be described below. All of 
these, with only two exceptions, lie outside the 
region containing the successful crystals or 
within 1° of the nearest boundary. Hence, a 
range of successful growth certainly exists. 
Single crystals of zinc plus 0.26 percent Cd can 
be grown inside it with complete success while 
outside only failures result. 

In Table I is given a classification of the 
specimens and the number of each type obtained. 
The successful crystals are placed in one main 
group A and the failures in three as follows: B, 
failure by change from single crystal of one 
orientation to single crystal of another orienta- 
tion; C, failure by being mosaic with no change 
of orientation ; D, failure by change from mosaic 
of one orientation to mosaic of another. The 


1 Because the crystals are very soft it is not easy to 
cleave them without some distortion. This makes it very 
difficult to determine the orientation closer than within 1°. 
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subdivisions of these groups are based on the 
character of the specimen (1) at the junction 
between nucleus and new crystal, (2) in the 
following section two or three cm in length, (3) 
in the rest of the specimen. Thus type Aa is 
single at the junction and remains so throughout ; 
type Ab is mosaic at the junction, then becomes 
and remains single; Bb is mosaic in the junction, 
becomes definitely single and then changes 
“‘discontinuously”’; etc. For a crystal to be 
mosaic at the junction means that the part of 
the original nucleus (always a single crystal) 
which is melted and then resolidified turns to 
mosaic (of the same orientation). Thus the new 
crystal gets its start as a mosaic.'® 

The types of change which may occur during 
growth are as follows: (1) “discontinuous,” a 
new crystal of different orientation starts and 
gradually crowds out the original orientation, the 
specimen growing as a double crystal for several 
cm; this was called a “sudden” change by 
Hoyem and Tyndall; (2) “slowly continuous,” 
all cleavage planes within the region of change 
(at least several cm) are apparently perfect, the 
orientation changing progressively by imper- 
ceptible amounts without the invasion of the 
crystal by a foreign orientation, the ‘‘gradual”’ 
change of Hoyem and Tyndall; (3) ‘‘rapidly 
continuous,” apparently the same as the pre- 
ceding type but occurring within 1 or 2 mm; 
(4) “typical evolution,” a type of mosaic change 
previously described by the writer.” 

Of the 75 group A crystals 68 were either 
single from the very start or became single very 
soon thereafter. Both these types (Aa and Ab) 
were distributed in the range quite at random. 
There were seven specimens (Ac) which remained 
mosaic for several cm below the start becoming 
gradually less mosaic (i.e., showing fewer and 
larger mosaic areas), finally becoming and re- 
maining single thereafter. These were all grown 
very near the range boundary. 

The appearance of the mosaic in the junction 
seems not to be connected in any way with the 
range, the percentages of mosaic junctions inside 
and outside being about equal. The phenomenon 
probably depends upon the precise conditions 


‘® Compare with work of Cinnamon and Poppy, refer- 
ence 3. 
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nucleus and the actual beginning of growth. 

Of the 129 failures (classes B, C, D) 55 started 
as mosaics but only 3 subsequently became 
single. Thus while all of the mosaic starts became 
single inside the range, practically none did 
outside. To start mosaic outside the range means 
to stay mosaic. Moreover, of the 74 failures 
which started single outside the range 27 (36 
percent) became mosaic thereafter. 

Fig. 1 shows that outside the range the region 
6>50° contains almost exclusively only group B 
failures, and that the region 6 <50° has relatively 
few. In the latter region group C and D specimens 
seem to be distributed fairly uniformly. 

The directions and magnitudes of the changes 
of orientation are also worth noticing. For @>50° 
all changes are either toward lower orientations, 
or else to 0@=90° or very near it. For 6<50° all 
changes are toward higher values of 6 and are 
either small (less than about 20°), or else all the 
way to 90° (i.e., for discontinuous changes only). 

It will be recognized at once that the writer's 
range resembles the H—T range in general char- 
acter though considerably distorted. Both ranges 
are narrow, have upper and lower temperature 
limits with maxima for an orientation near 35° 
and extend upward indefinitely for @ near 90°. 
While both ranges determine regions of successful 
growth, the types of change of orientation or 
failures are more varied and complex in the 
present work.'® In two respects, however, the 
present results are very different. First, there are 
three regions or sets of orientations (in the 
neighborhood of 22°, 45° and 90°) in which 
successful crystals can be grown over a very 
wide range of temperatures. Second, it appears 
impossible to grow crystals having orientations 
between about 50° and 85°, whereas one of the 


16 It was early suspected that perhaps H-T had also 
grown mosaics but had not detected them and had counted 
them as singles. As examples of their ‘“‘failures’’ some 15 
short sections of specimens in which a sudden change had 
occurred were preserved. Recent examination of these 
showed that 7 were of the type Ba while the other 8 were 
good singles changing to mosaic just before a discontinuous 
change, i.e., these latter do not belong exactly to any of 
the writer's classes. An insufficient number of successful 
crystals were preserved to warrant the drawing of any 
definite conclusion. However, of nine such crystals or 
pieces of crystals, only two showed any indication of 
optically mosaic structure and both of these were grown 
near the border of the range. 
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most outstanding features of the H-T range is 
the wide latitude of growth conditions for single 
crystals having orientations above about 65°. 
Experiments have also been done at the same 
speed of growth (1.2 cm/min.) with other 
concentrations (0.00 percent, 0.08 percent, 0.20 
percent, 38 percent, 53 percent, 0.99 percent Cd) 
though the whole region of orientations and 
temperatures was not completely explored in 
each case. Throughout this work only truly single 
crystals containing at least 0.20 percent Cd were 
used as nuclei.’ For pure zinc only type Cb 
specimens were obtained (mosaics with no per- 
ceptible change of orientation) over widely 
different temperature conditions. Thus we have 
for pure zinc a sort of stability of growth of 
mosaics. After adding 0.08 percent Cd, a few 
single crystals were obtained at low temperatures 
(about t—t=13°). All other specimens were of 
type Dc or Dd. There was no indication of a 
range. The results for 0.20 percent and 0.38 
percent Cd were in the main the same as for 
0.26 percent Cd, though the ranges obtained 
were narrower. For instance, in the case of 0.38 
percent the range was only about half as wide 
at 6=42° as for 0.26 percent Cd. Apparently 
the Cd concentration yielding the maximum 
width of the range is somewhere near 0.26 
percent. For 0.53 percent Cd the range was 
confined to a narrow strip near 90° orientation 
since crystals successfully started at other ori- 
entations changed soon (type Ba) to singles of 
nearly 90° orientation. It was found furthermore 
that with increase of Cd content (from 0.08 
percent up) mosaics and mosaic changes became 
less frequent and never occurred for 0.99 percent 
Cd plus Zn. Continuous and discontinuous 
changes from single crystal of one orientation to 
single of another become, on the other hand, 
more frequent (outside the range, of course). 
With 0.99 percent Cd singles could not even be 
started.'§ Several crystals of different orienta- 


1 This was necessary since specimens of Cd pure zinc 
were mosaic and unsuitable for use as nuclei. It was also 
hoped by this use of nuclei containing Cd to start single 
crystals even when using pure zinc and thus discover a 
range of growth for it. This device, though successful with 
another method of growth (Poppy, reference 3) failed for 
the C-G method. 

8 Cd is soluble in Zn to the extent of 1 percent at room 
“<n Circular of Bureau of Standards, No. 395, 
p. 140. 
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Fic. 2. Range obtained with 0.26 percent Cd+Zn at 
— of 0.55 cm/min. Symbols have the same meaning as 
in Fig. 1. 


tions always started at the junction. These went 
through all sorts of changes of orientations, both 
continuous and discontinuous, to higher and 
lower orientations, rapidly and slowly. 

These findings for 0.99 percent Cd in Zn seem 
to agree with those of Rosbaud and Schmid!® 
who also found it impossible to grow single 
crystals of Zn+Cd by the C-G method when 
the Cd concentration was above 1.03 percent. 
Straumanis,”° however, obtained specimens con- 
taining as much as 3 percent Cd which had all 
the appearance of single crystals and which he 
therefore called single. These were grown by the 
earlier method of Bridgman”! and at the very 
low speed of 1.3 cm/hr. 

The 0.26 percent Cd+Zn was studied also at 
the speed of 0.55 cm/min. The results are shown 
in Fig. 2. While the field was not explored as 
thoroughly as for the higher speed, it is clear 
that a range certainly exists although not the 
same range as before. The two most noticeable 
differences are first a shift of about 5° of the 
left-hand branch of the range and a widening 
out of the center branch sufficiently to join it 
definitely to the strip for 6 near 90°. Lowering 
the speed, therefore, increases tremendously the 
possibility of growth of single crystals of this 


19 P, Rosbaud and E. Schmid, Zeits. f. Physik 32, 197 
(1925). They do not state what speed of growth was used. 

20 M. Straumanis, Zeits. f. anorg. allgem. Chemie 180, 1 
(1929). 

21 P, W. Bridgman, Proc. Am. Acad. 60, 308 (1925). 
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combination of Zn and Cd. These results are 
different from those of Hoyem and Tyndall who 
found that increasing the speed definitely raised 
the lower limit of their range. Nor are they in 
agreement with the hypothesis of Cinnamon” to 
the same effect, which, however, may not apply 
to impure zinc nor to growth of much smaller 
diameter crystals by another method at a 
different speed. The changes of orientation are 
similar to those for the higher speed, except that 
they are less in magnitude, rarely exceeding 7° 
or 8°. 

Rosbaud and Schmid" and Straumanis,”’ in 
growing crystals of zinc with Cd impurity, found 
that the Cd entered the growing crystal in the 
same proportion in which it was present in the 
melt. This was found to be the case in the 
present work also.** They found, however, that 
the Cd did not distribute itself uniformly in the 
crystals but was segregated in layers parallel to 
the basal plane of the hexagonal lattice. This 
was evidenced in Straumanis’s experiments by 
the appearance of uniformly spaced ridges on 
the surfaces of his crystals parallel to that plane. 
These appeared for Zn containing Cd (from mere 
traces up to 5 percent Cd**) but not redistilled 
and therefore presumably Cd free Zn. Etching 
experiments led Straumanis to the conclusion 
that the ridges were indeed richer in Cd than 
the layers between them. 

The present writer finds the same periodic 
ridges*®> for Cd concentrations between 0.20 


“= Cinnamon’s hypothesis* states that the “ratio of the 
temperature gradient to the rate of growth must be near 
a certain favorable, but not necessarily the same, value for 
each orientation.’ According to this one should expect to 
find for lower speeds a range of successful growth at a lower 
temperature gradient, and therefore at a lower temperature 
in the present work in which f—¢) is proportional to the 
gradient. 

*3 By spectroscopic analysis carried out by Professor 
Tyndall. 

** Straumanis’s specimen containing 5 percent Cd was a 
polverystal. However, each component showed the same 
typical periodicities in Cd segregation as did his ‘‘single”’ 
crystals containing 3 percent Cd or less. 

** The fine drawings and photographs by Straumanis, 
some of which were reproduced recently in a paper by 
F. Zwicky, Rev. Mod. Phys. 6, 193 (1934), illustrate 
completely and in detail the appearance of the present 
writer's crystals. 
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percent and 0.99 percent. They were prominent 
enough to be easily visible to the naked eye, and 
their prominence increased with the Cd content. 
Because of this and because they were always 
exactly parallel to the primary cleavage planes 
they made possible the determination of the 
orientation of the crystals before cleavage. By 
means of them it was quite easy also to follow 
the course of discontinuous and slowly continu- 
ous changes of orientation merely by examination 
of the crystal surfaces. 

It seems quite certain in the light of the 
foregoing that the existence of the H—-T range 
must be attributed mainly to the presence of 
the Cd. That the present work has not exactly 
duplicated it is probably due to the fact that 
other impurities besides Cd were present in the 
H-—T zine and undoubtedly exerted a modifying 
influence. It would be interesting to test this 
point directly by adding other elements besides 
the Cd in known amounts. This has been tried 
in a preliminary way only for Pb. The addition 
of only 0.02 percent Pb to 0.26 percent Cd+Zn 
apparently leaves the range unchanged, while 
0.2 percent Pb very materially modifies it. 

There is nothing in the above to disprove the 
hypothesis that crystals of pure (99.99 percent) 
zinc must have certain favorable growth condi- 
tions which are a function of the orientation. 
These conditions simply cannot be investigated 
by the C—G method?* on account of the appear- 
ance of mosaics. 

This work was done at the suggestion and 
under the direction of Professor E. P. T. Tyndall. 
For his very great help, constant interest and 
encouragement I thank him most heartily. 

Most of the work was done in the writer's 
laboratory at Union College, Lincoln, Nebraska. 
Apparatus was loaned by the departments of 
Physics of the State University of Iowa and the 
University of Nebraska. I am also very grateful 
for this help. 


26 That is, with any technique which the writer has tried 
so far. Given a way to eliminate mosaics there is no reason 
why the C-G method should not be used as a means of 
studying growth of single crystals of pure zinc. 
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In previous papers we have investigated the conditions for spontaneous division of systems, 
which are seats of physico-chemical reactions of any kind. In this paper the previous general 
results are specialized for one particular set of reactions: oxidation of glucose and glycolysis. 
It is found, that whenever the two occur simultaneously, the ratio of their relative rates de- 
termines whether the system will divide spontaneously or not. Increased glycolysis favors 
spontaneous division. It is pointed out that the rapidly dividing cancer cells are characterized 


by an abnormally high rate of glycolysis. 





N a paper, published several years ago in this 

journal' we have shown how various general 
phenomena in living cells are in a natural way 
brought within the scope of studies of mathe- 
matical physicists. Such studies open a large 
field for physico-mathematical investigations, 
confronting us with problems, which being 
essentially physical in their nature, have never- 
theless been very little, if at all, considered 
hitherto by physicists. 

In subsequent papers*~° we have investigated 
mathematically various aspects of spontaneous 
division of small systems, which are seats of 
physico-chemical reactions. In the future we 
shall for brevity call such systems ‘‘metabolizing 
systems.’ Typical representatives of such sys- 
tems are living cells. In fact metabolism is the 
only physical characteristic common to all grow- 
ing cells,*:* though they are of course more 
complex than the systems which we hitherto 
studied. 

In calculating various factors which contribute 
to spontaneous division, we considered only such 
cases in which only one substance is being pro- 
duced or absorbed in the system. When the 
system metabolizes a great number of different 
substances, the resulting effects are given by the 
sum of the contributions of all individual sub- 
stances. The condition that above a critical size 
the system should become unstable and eventu- 
ally divide is now expressed by 


‘ Rashevsky, Physics 1, 143 (1931). 

* Rashevsky, Physics 2, 303 (1932). 

* Rashevsky, Physics 5, 374 (1934). 

* Rashevsky, Physics 6, 33 (1935). 

* Rashevsky, Physics 6, 35 (1935). 

® Philosophy of Science 1, 176 (1934). 


* Cold Spring Harbor Symposia on Quantitative Biology 
2, 188 (1934). 


Ygila(RT/M—T’)/D.—(T—-1")/Di}k>0, (1) 
k 


the sum being taken over all the metabolized 
substances. Eq. (1) replaces the simpler Eq. (47) 
of the former paper,’ the notations being the 
same. 

In the present paper we shall investigate a 
specific case of a set of reactions, frequently 
found in nature. ‘ 

If we have a system in which all the reactions 
can be separated in two groups, one characterized 
by relatively large g,’s and the other by much 
smaller ones, we may in the first approximation 
extend the summation only over the first group. 
Such a situation actually is found in most living 
cells. For reactions connected with respiration 
qg~ 10-4 g-cm~*- min... Those reactions include : 
Consumption of oxygen (Oz) and glucose (g); 
production of carbon dioxide (CO:), lactic acid 
(1) and water (w). Considering the next im- 
portant reactions, those connected with protein 
metabolism, we find g~10-7 g-cm-*-min.~. 
Others are still slower. 

It is of special interest therefore to calculate 
the resulting effect of the five above-mentioned 
reactions. We do this by extending in (1) the 
summation over indices: ;, 2, 3, 4, 5, with the fol- 
lowing meaning: 


gi=qi>0; ge=gceo2>0; gs=qu>0; 


Gs=9 <0; Gs=qo2<0. (2) 


Indices of the other constants of (1) have 

correspondingly similar meanings. 
Multiplication of all D, and D; with the same 

factor does not change (1); the same holds about 
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the g’s. Hence we need to know only the relative 
values of the rates of reactions as well as of the 
diffusion coefficients. 

Since the diffusion coefficient of the dissolved 
substance is determined by the resistance offered 
by the viscosity of the solute, and since approxi- 
mately the latter may be calculated by the 
application of Stokes’ formula,* considering the 
dissolved molecule as spherical, D must vary 
approximately as the inverse of the third root 
of molecular volume, which latter in its turn 
varies approximately as the molecular weight. 
We could of course introduce exact atomic 
volumes, but it would only somewhat complicate 
the calculations without changing the conclu- 
sions. The use of Stokes’ formula is only roughly 
approximate, and does not justify greater exacti- 
tude in this next step. Substituting D~1/M! 
into (1), and expressing, for reasons which will 
become obvious immediately, the relative rates 
of reactions not in g-cm~*-min.~', but in mole- 
cules* per cm* per min., we find 


Yailaa(1—1’M/RT) M3 


—b(T—1’)M*3/RT},>0, (3) 


a and b being coefficients of proportionality in the 
formulae, connecting D and 1/M!. 

The five reactions considered here are con- 
nected together in the following way. Some 
molecules of glucose are simply split (glycolyzed) 
into lactic acid, according to the equation: 


CgH 204 = 2C3H Os. (4) 


Some are oxidized to CO, and water, according 
to 


CH 20. + 602 = 6CO2+6H,0. (5) 


It is customary to characterize the metabolism 
of a cell by the ratio 8 of the two reactions.’ 
Then 38 represents the ratio of the number of 
glycolyzed molecules of sugar to the number 
completely oxidized.* If one molecule of glucose 
is oxidized, then 38 are glycolyzed. Altogether 
38+1 glucose molecules are consumed. For the 


SE. Baars, Geiger-Scheel Handb. d. 
(Springer, Berlin, 1928). 

* That is substituting for g, 9M) N, N—being Avogadro's 
number. 

*O. Warburg, The Metabolism of Tumors, Constable 
and Company, London, 1930. 
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oxidation of one molecule 6 molecules of Os are 
necessary. 68 molecules of lactic acid and 6 
molecules of COs are produced. 

Hence the relative values of the g’s are: 


qi = 68 ; q2=6; 93 =6; gs= — (3841); q5= —6. (6) 


An approximate estimation has shown us’ 
that T and I” may be of the same order of 
magnitude as R7/M. A more exact calculation 
which will be given elsewhere, shows that 


r=L/2M, (7) 


where L denotes the work necessary to bring 
one gram mole of dissolved substance from an 
ideal gas state into solution. A similar expression 


holds for I’. 
Substituting (7) into (3) we find: 
Yai aa(t—L'/2RT)M" 
1 
—b(L—L’)M*?/2RT],>0, (8) 








or, denoting the expression in square brackets by 
A; and substituting (6) into (8): 


(6A \8+6A2+6A3)/(3A,s8+A4+6A;)>1, (9) 


which expresses the condition of spontaneous 
division of the system in terms of 8. 

For solutions of oxygen and COs, in water the 
L’s can be calculated from the Bunsen partition 
coefficients and are found to be: 


Leog= —1-7X10" ergs; 
(10) 
Lo2= —3-8 X10" ergs, 


the minus sign indicating a force of repulsion 
between the molecules of the solvent and solute.’ 
For the determination of L’ no data are available 
at all. Nor do we know anything about L for 
glucose and lactic acid. Of the a’s we know only 
that they are of the order of unity. 

Putting L,.=L,’, ag=1, which should hold at 
least as a first approximation, and introducing 
(6) and (10), together with the numerical values 
of M’s and RT into (8) we find: 


(278+ 45) /(178+40) >1. (11) 


The left-hand sides of (9) and (11) represent the 
ratio of the sum of the positive terms of (8), 
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respectively, (1), to the sum of the negative ones. 
But this is nothing else* than the ratio of the 
sum of ‘‘dividing’’ forces, which satisfy Eq. (47) 
of the previous paper,’ to the sum of the forces 
opposing division. From (11) we see that this 
ratio is the larger, the larger 8. It is interesting 
to note in this connection that O. Warburg® has 
found that the abnormally easily dividing cancer 
cells are characterized by an abnormally high 8. 
Whereas for cells of normal tissues 8 is practically 
zero, for cancer cells it is as high as 4. With 
8=0 the ratio (11) is 1-12; with 8B=4 it is 1-42. 

It may be questioned whether water must be 
included amongst the substances produced in 
the system, if the latter itself consists mainly of 
water, as in the case of living cells. If the latter 
grow at such a rate as to double in about 10-24 
hours, then all the water produced remains in 
the system to build it up and does not diffuse 
outwards. Such substances do not contribute to 
the osmotic and other forces.':? * Omitting the 
term corresponding to water in (8) and (9), that 
is putting g;=0, we find 


(278+30) /(178+40) >1, (12) 


which is satisfied only if 8>1. For 8=0 no 


division can occur. For 8=4, however, the 
dividing forces exceed the opposite ones by 28 
percent. 

If we consider only the effect of osmotic 
forces, that is if we put T,=I,,’=0, then (8) 
leads to 


(278+21)/(178+25) >1, (13) 
which shows that for 8=0 no division can occur 
at all. (13) is satisfied only for 8>0-36. But for 
8=4 the ratio (13) equals 1-39. 

The physical reason for this can easily be 
visualized by remembering* that produced sub- 
stances exert an osmotic pressure directed out- 
wards and thus tend to divide the system into 
smaller parts. Consumed substances produce 
osmotic forces in opposite directions. Strong 
glycolysis results in an increase of production of 
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lactic acid and in a relative decrease of con- 
sumption of oxygen. 

The situation will not be changed materially 
by taking into account I for glucose and lactic 
acid. Though no direct data to calculate it are 
available, an indirect estimation of the effect of 
such a correction can be made. Glucose in solu- 
tions slightly increases the surface tension of 
water. Hence there is a force of attraction 
between molecules of glucose and of water. I is 
positive and its effect is therefore opposed to 
that of osmotic forces. It will tend, therefore, to 
decrease the ‘‘restoring” forces due to glucose. 
For lactic acid no data on surface tension are 
available, but from analogy to other similar 
compounds and from general theoretical con- 
siderations we should expect it to be capillary- 
active. This means a negative I and enhances the 
dividing forces produced by lactic acid. In any 
case, however, relations like (12), (11) or (13) 
are obtained. 

We thus see that making any assumptions 
within a plausible range about various uncer- 
tainly known physical factors leads us to the 
same result, namely: Increased glycolysis leads 
to an increase of dividing forces and thus favors 
spontaneous division. 

When (8) or (9) are satisfied, the critical size 
at which a spontaneous division occurs can be 
calculated in a similar manner as before* by 
equating the left-hand side of (1) to 8y/R,?. 
The values thus obtained are of the order of 
magnitude of 10-* cm. 

The calculations in this paper were based on 
(1) which holds for single-phase systems. How- 
ever a similar expression is derived also from the 
formulae which we developed for two-phase 
systems,‘ and the general conclusions remain 
the same. 

We have also studied the electric charges 
which are produced by metabolism.* Their 
possible effects on spontaneous division will be 
discussed elsewhere. Except for some special 
cases they are found to be of secondary impor- 
tance as compared with those studied here. 
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